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ABSTRACT  OF  THE  DISSERTATION 


Shelf  Circulation  in  the  Gulf  of  California 

by 

Mark  A.  Merrifield 

Doctor  of  Philosophy  in  Oceanography 

University  of  California,  San  Diego,  1989 

Professor  Clinton  D.  Winant,  Co-chairman 
Professor  Robert  T.  Guza,  Co-chairman 


The  first  long-term  moored  observations  in  the  Gulf  of  California  (1983-84)  are  used  to 
describe  the  dynamically  important  spatial  and  temporal  scales  of  the  shelf  circulation  on  oppo¬ 
site  sides  of  the  150-km  wide  gulf,  and  to  investigate  the  characteristics  of  coastal-trapped  waves 
that  propagate  into  the  gulf  along  the  mainland  shelf.  A  description  of  the  observations  is  con¬ 
tained  in  chapter  1.  The  main  results  of  this  section  are: 

1)  Local  wind  forcing  is  not  observed. 

2)  Flow  on  the  two  facing  shelf  regions  is  uncorrelated,  although  recirculation  with  basin  flow 
is  observed  on  both  shelves.  The  net  surface  transport  is  correlated  with  the  the  cross-gulf  pres¬ 
sure  gradient. 

3)  A  rise  in  subinertial  sea  level  on  the  mainland  shelf  is  accompanied  by  a  drop  in  isotherm 
depths  such  that  pressure  fluctuations  rapidly  decay  with  depth. 

4)  Cross-shelf  geostrophy  is  weakly  apparent  on  both  shelves.  The  acceleration  of  alongshclf 
currents  is  correlated  with  the  alongshelf  pressure  gradient  on  the  mainland  shelf,  indicating  the 
presence  of  remotely  forced  waves.  Only  temperature  observations  show  evidence  of  waves 
reaching  the  Baja  California  peninsula. 

In  chapter  2,  an  empirically  derived  propagating  signal  is  compared  with  a  coastal -trapped 
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wave  model.  The  model  is  forced  by  a  simple  extrapolation  of  storm  winds  to  the  coast.  The 
main  results  of  this  chapter  are: 

1)  The  observed  temporal  behavior  of  wave  events  is  consistent  with  storm  forcing  of  mode  1 
coastal-trapped  waves. 

2)  The  observed  and  mode  1  density  fields  are  similar  in  structure,  although  observed  ampli¬ 
tudes  are  nearly  three  times  larger  than  predicted.  The  observed  offshore  pressure  field  decays 
more  rapidly  with  depth  and  has  a  shorter  offshore  decay  scale  than  mode  1. 

3)  Alongshelf  currents  are  barotropic  during  energetic  wave  events  at  the  midshclf,  in  accor¬ 
dance  with  mode  1.  Mode  1  cannot  account  for  the  observed  cross-shelf  structure  in  the 
alongshelf  flow  field. 

4)  Higher  modes  are  not  observed. 

In  chapter  3,  complex  empiric  J  orthogonal  functions  arc  assessed  as  a  possible  method  for 
isolating  propagating  variability  in  the  gulf  data  set.  The  results  of  this  study  indicate  that 
CEOFs  are  inadequate  for  the  detection  of  the  wave  signal  due  to  the  wide-banded  frequency 
characteristics  of  the  signal. 
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CHAPTER  1 


SHELF  CIRCULATION  IN  THE  GULF  OF  CALIFORNIA: 
A  DESCRIPTION  OF  THE  VARIABILITY 


1.  INTRODUCTiON 

Recent  studies  of  shelf  circulation  have  identified  important  dynamic  processes  of  open 
coasts,  e.g.,  cross-shelf  geostrophy,  upwelling,  and  coastal-trapped  waves  (see  Brink,  1987  for  a 
review).  Circulation  along  sheltered  coasts,  as  in  a  marginal  sea,  has  received  less  attention. 
Marginal  seas  offer  a  particular  advantage  for  shelf  studies  in  that  open  ocean  influences  are 
diminished  or  more  easily  monitored.  In  turn,  an  understanding  of  shelf  circulation  in  a  marginal 
sea  may  give  insight  about  the  circulation  of  the  basin  as  a  whole  since  the  shelves  constitute  a 
significant  percentage  of  the  total  basin  area.  The  Gulf  of  California  is  similar  to  a  number  of 
other  marginal  seas  (Adriatic  Sea,  Red  Sea,  Persian  Gulf)  in  that  it  is  long  and  narrow  (1000  km  x 
150  km),  and  open  at  the  southern  boundary.  From  1983-1984,  the  gulf  was  the  subject  of  a  joint 
field  experiment  between  Scripps  Institution  of  Oceanography  (SIO)  and  Centro  de  Investigation 
Cicntlfica  y  de  Education  Superior  de  Ensenada  (CICESE).  The  observations  in  the  Guaymas 
Basin  (Fig.  1),  are  the  first  long-term,  moored  measurements  of  currents,  temperature,  and  bottom 
pressure,  as  well  as  winds,  made  in  the  Gulf  of  California.  The  objectives  here  are  to  present  a 
basic  description  of  these  observations  on  seasonal,  low  frequency,  and  tidal  time  scales,  to  com¬ 
pare  and  constrast  the  two  facing  shelf  regions,  and  to  provide  a  qualitative  assessment  of  the 
relative  importance  of  local  and  remote  forcing  on  gulf  shelf  circulation. 

Of  particular  relevance  to  this  study  are  the  sub-tidal  sea  level  events,  with  positive 
elevations  of  10-30  cm  and  time  scales  of  about  10  days,  which  propagate  poleward  along  the 
Mexican  coast  into  the  Gulf  of  California  (Christensen  et  al.,  1983,  Enfield  and  Allen,  1983). 
Enfield  and  Allen  showed  that  the.co  events  are  associated  with  tropical  storms  and  hurricanes  that 
occur  during  the  summer-fall  season  off  southern  Mexico.  Phase  speeds  indicate  forced  propaga¬ 
tion  in  the  region  of  storm  activity  slowing  to  free  propagation  in  the  gulf,  where  they  are 
somehow  dissipated,  as  they  aie  not  detected  on  the  Pacific  Ocean  coast  past  the  Baja  California 
peninsula.  Simple  estimates  of  phase  speed,  length  scale,  and  mass  transport,  based  exclusively 
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on  sea  level  observations,  suggest  that  these  events  are  hybrid  form  coastal-trapped  waves,  more 
similar  to  internal  Kelvin  than  continental  shelf  waves.  The  present  study  provides  information 
about  the  temperature  and  current  fluctuations  associated  with  these  waves. 

A  brief  review  of  circulation  in  the  Gulf  of  California  is  given  in  section  2,  and  the 
experiment  is  described  in  section  3.  In  section  4,  the  basic  observations  (wind,  temperature, 
pressure,  and  currents)  are  presented,  with  emphasis  on  identification  of  spatially  coherent  pat¬ 
terns.  Rough  momentum  balances  are  given  in  section  5,  a  coastal-trapped  wave  event  is 
described  in  section  6,  and  a  summary  follows  in  Section  7. 


2.  BACKGROUND 

The  Gulf  of  California,  or  Sea  of  Cortes,  is  located  on  the  west  coast  of  Mexico  between 
the  Baja  California  peninsula  and  the  Mexican  mainland  (Figure  1).  The  gulf  is  divided  into 
three  distinct  topographic  regions.  The  northern  region  is  shallow  (typically  <  200  m)  and 
comprises  roughly  a  third  of  the  gulf’s  length.  The  southern  province,  in  open  communication 
with  the  Pacific  Ocean  to  depths  of  3000  m,  is  trough-shaped  with  a  series  of  deep  basins  (1000- 
3000  m)  flanked  by  a  steep,  rocky  Baja  California  shelf  and  a  wider  mainland  shelf.  Numerous 
islands,  capes,  and  inlets  arc  found  on  either  shelf.  A  third  region  of  midriff  islands,  narrow  pas¬ 
sages,  and  shallow  sills  separate  the  northern  and  southern  gulf. 

Atmospheric  forcing  over  the  gulf  is  characterized  by  weak  summer  winds  towards  the 
northwest  and  stronger  winter  winds  towards  the  southeast  (Roden,  1964).  Winds  are  polarized 
along  the  gulf  axis  by  an  uninterrupted  mountain  range  along  the  Baja  California  peninsula 
(elevations  higher  than  700-1000  m)  and  to  a  lesser  extent  by  the  Sierra  Madre  Occidental  (higher 
than  1500  m)  further  inland  on  the  Mexican  mainland.  The  peninsular  mountains  effectively  iso¬ 
late  the  gulf  from  the  moderating  influence  of  the  Pacific  Ocean  resulting  in  a  hot,  arid  climate 
more  typical  of  the  surrounding  desert  regions  than  the  Pacific  coastal  zones.  Conditions  are 
somewhat  more  tropical  during  the  summer  when  prevailing  winds  advect  moist  air  in  from  the 
south.  Tropical  storms  and  hurricanes  are  active  from  May  to  October  in  the  eastern  north 
Pacific.  These  energetic  disturbances,  with  wind  speeds  exceeding  100  knots,  form  at  low  lati¬ 
tudes  and  travel  northwestward  along  the  Mexican  Pacific  coast.  The  majority  of  storms  turn 
westward  before  reaching  the  gulf,  a  lesser  number  turn  toward  the  coast,  while  very  few  enter 
the  gulf  region  directly  (Roden.  1964).  Only  one  storm  crossed  the  southern  mouth  of  the  gulf 
during  this  experiment. 
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The  circulation  of  surface  waters  in  the  southern  gulf  was  first  described  by  Thorade 
(1909),  using  ship  logs,  as  cyclonic  with  net  surface  transport  in  the  direction  of  seasonal  winds. 
He  also  related  colder  mainland  coast  temperatures  during  winter  to  upwelling  favorable  winds. 
The  cyclonic  circulation  pattern  was  also  seen  in  late  summer  in  hydrographic  data  (Roden  and 
Groves,  1959)  with  low  salinity  water  flowing  northward  on  the  mainland  coast  and  high  salinity 
water  flowing  southward  along  the  peninsula  at  50-m  depth.  A  narrow  jet  of  high  salinity  outflow 
was  also  seen  along  the  Baja  California  coast  at  the  gulf  entrance  during  December  (Roden, 
1972)  suggesting  that  the  outflow  may  extend  along  the  peninsula.  Bray  (1988),  in  a  study  of 
thermohaline  circulation  in  the  gulf,  analyzed  geostrophic  currents  at  the  same  cross-gulf  transect 
as  our  moored  array.  Her  findings  show  cyclonic  circulation  in  the  Guaymas  Basin  during  sum¬ 
mer,  weaker  anti-cyclonic  flow  during  spring  and  fall,  and  the  vertical  structure  of  the  flow  aver¬ 
aged  across  the  gulf  consisting  of  outflow  between  50  m  and  250  m,  inflow  between  250  m  and 
500  m,  and  a  50-m  deep  surface  layer  where  transport  is  in  the  direction  of  seasonal  surface 
winds. 

The  above  description  is  representative  of  the  basin  scale  flow.  On  shorter  scales,  circu¬ 
lation  may  be  more  complex.  Brown  (1965),  in  a  series  of  drogue  studies  in  the  southern  gulf, 
reported  an  irregular  flow  field  changing  with  geographic  location  and  depth.  Eddy-like  features 
have  been  observed  in  the  Guaymas  Basin  with  drogues  (Emilsson  and  Alatorre,  1980)  and  in  sea 
surface  temperature  (e.g.,  Badan-Dangon  et  al,  1985).  Events  which  resemble  the  squirts 
observed  off  northern  California  (Kelly,  1985)  have  been  detected  north  of  Guaymas  using  dro¬ 
gues,  current  meters,  and  satellite  images  (M.  Lavin,  CICESE,  unpublished  data,  1988).  Cold  sea 
surf?  ;e  temperature  plumes  have  been  reported  on  both  sides  of  the  gulf  and  related  to  the  onset 
of  upwelling  favorable  winds  (Badan-Dangon  et  al.,  1985). 

Tides  in  the  gulf  are  mixed  with  considerable  difference  between  the  diurnal  and  semi¬ 
diurnal  components.  Filloux  (1973),  using  observations  of  tidal  elevation  on  the  perimeter  of  the 
gulf,  described  the  semi-diurnal  tide  as  progressive  towards  the  north  with  amplitudes  in  the  cen¬ 
tral  gulf  of  0.1 -0.3  m,  increasing  rapidly  near  the  midriff  islands  to  nearly  2  m  at  the  head  of  the 
gulf.  In  contrast,  the  diurnal  tide  behaves  like  a  standing  wave  with  amplitudes  ranging  from 
0.25  at  the  entrance  to  0.5  m  at  the  head  of  the  gulf.  The  large  semi-diurnal  tidal  range  in  the 
northern  gulf  is  accompanied  by  strong  tidal  currents,  particularly  in  the  island  passages  where 
speeds  of  3  m/s  have  been  observed  (Roden  1964).  Hendershott  and  Speranza  (1971)  proposed 
that  the  strong  currents  in  the  north  give  rise  to  large  tidal  dissipation,  manifested  as  a  virtual 
amphidromc  for  the  M2  tide  in  the  central  gulf.  The  strong  tidal  flows  also  generate  internal 
waves  which  contribute  an  estimated  10%  to  the  dissipation  of  the  semi-diurnal  tide  (Fu  and  Holt, 
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1984). 

3.  EXPERIMENT 

The  principal  study  site  was  a  transect  between  Guaymas  on  the  mainland  coast  and 
Santa  Rosalia  on  the  Baja  California  peninsula,  in  the  guaymas  Basin  of  the  central  gulf  (Figure 

1) .  The  instruments  were  deployed  in  four  phases,  each  approximately  six  months  long  (Figure 

2) .  In  the  first  phase  (11/82-5/83),  currents  were  measured  on  the  Guaymas  shelf  to  determine 
typical  current  speeds  in  this  region  prior  to  more  extensive  deployments.  The  second  phase 
(5/83-11/83)  was  a  pilot  experiment  designed  to  compare  currents  separated  in  the  cross-gulf 
(Guaymas-Santa  Rosalia)  and  along-gulf  (Guaymas-Topolobampo)  directions.  Bottom  pressure 
measurements  were  also  made  at  these  three  shelf  locations.  The  spatial  scales  of  lower  atmos¬ 
pheric  variables,  such  as  winds,  pressure,  and  temperature,  were  sampled  during  July  1983  by 
extensive  aircraft  overflights  in  the  northern  half  of  the  gulf  (Candela  et  al.,  1984).  During  the 
third  phase  (11/83-5/84),  fixed  PAMII  (Portable  Automated  Mesonet)  weather  stations  were 
deployed  to  measure  winds,  air  temperature,  and  barometric  pressure  continuously  for  the 
remainder  of  the  study  (Merrifield  et  al.,  1986).  A  second  set  of  overflights  was  made  during 
March  1984  (Candela  et  al.,  1985).  Moorings  were  maintained  at  Guaymas  and  Santa  Rosalia  to 
observe  winter  conditions  on  the  two  shelves.  The  main  deployment  (5/84-11/84)  included  a 
denser  array  across  the  Guaymas-Santa  Rosalia  transect  with  an  additional  mooring  in  the  deep 
basin  near  the  Guaymas  slope,  a  mooring  at  Topolobampo,  and  pressure  sensors  at  all  three  shelf 
locations  as  well  as  at  San  Francisquito  Bay  and  Isla  Tiburon,  located  north  of  the  transect.  In 
addition,  sea  level  observations  at  Salina  Cruz,  Acapulco,  Guaymas,  Santa  Rosalia,  and  Cabo  San 
Lucas  (Figure  1)  were  obtained  from  pre-existing  tide  gauges  maintained  by  UNAM  and 
CICESE.  Mooring  locations  are  shown  in  Figure  1;  the  depths  and  deployment  times  of  each 
sensor  are  listed  in  Table  1. 

Currents  were  usually  measured  (Table  1)  with  Vector  Measuring  Current  Meters 
(VMCMs,  Weller  and  Davis,  1980),  but  also  with  General  Oceanics  inclinometers  (Barton  et  al., 
1980).  Records  from  a  VMCM  and  inclinometer  at  70-m  and  75-m  depth  below  the  surface  on 
the  Guaymas  shelf  (M7)  were  in  only  fair  agreement  at  low  frequencies  (correlation  =  0.7)  with 
the  inclinometer  yielding  a  noiser  record.  Temperature  was  measured  at  each  current  meter  and 
pressure  sensor  package  and,  during  the  main  deployment,  with  thermistor  chains  in  200-m  depth 
on  the  Guaymas  and  Santa  Rosalia  shelves.  Bottom  pressure  was  measured  by  Paroscientific 
pressure  sensors  in  a  package  designed  and  fabricated  at  SIO  (Erdman,  1983).  A  complete 
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presentation  of  the  moored  observations  and  details  of  the  data  processing  are  given  in  Merrifield 
et  al.  (1986). 

Low  frequency  or  subinertial  time  series  in  this  paper  refer  to  series  which  have  been 
convolved  with  a  Groves  51-weight  filter  (Groves,  1955)  and  decimated  to  6-hourly  values.  The 
Groves  filter  has  a  half-power  point  at  approximately  0.3  cpd  and  is  chosen  for  its  ability  to 
suppress  leakage  from  strong  tidal  harmonics  and  its  relative  short  length  in  the  time  domain, 
which  minimizes  data  loss  in  time  series  with  long  gaps.  In  the  following,  a  significant  correla¬ 
tion  between  two  time  series  excedes  the  95%  significance  level.  For  the  typical  number  of 
degrees  of  freedom  (Davis,  1976),  this  corresponds  to  a  correlation  of  at  least  0.25. 


4.  DESCRIPTION  OF  OBSERVATIONS 

Winds 

Time  series  of  low  frequency  winds  at  each  of  the  PAMII  weather  stations  (Figure  3) 
illustrate  the  dominant  seasonal  pattern  of  winds  over  the  northern  half  of  the  gulf  with  winds 
toward  the  northwest  during  the  summer  and  toward  the  southeast  during  the  remainder  of  the 
year.  Winds  are  strongest  during  winter  with  maximum  speeds  at  Isla  Tortuga  (W5)  of  15  m/s 
compared  to  10  m/s  during  the  summer.  Principal  axes  of  the  low  frequency  winds  are  directed 
primarily  along  the  gulf  with  local  variations  in  orientation  probably  due  to  topographic  steering 
(Figure  4).  Standard  deviations  are  typically  2-5  times  larger  in  the  along-gulf  than  cross-gulf 
wind  component  The  exception  is  at  Isla  Piojo  (W2)  where  cross-gulf  winds  are  as  strong  as 
along-gulf  winds,  likely  due  to  lee  waves  over  the  Baja  California  mountains  (C.  E.  Dorman,  per¬ 
sonal  communication). 

The  spatial  coherence  of  the  low  frequency  wind  field  is  examined  with  empirical 
orthogonal  functions  (EOFs).  In  this  analysis,  U  and  V  wind  components  are  combined  into  a 
vector  time  series  and  EOFs  obtained  from  the  complex  covariance  matrix  (e.g.,  Kundu  and 
Allen,  1977).  This  method  is  preferable  to  a  scalar  analysis  because  wind  vectors  sometimes 
rotate  in  time,  with  the  U  and  V  wind  time  series  out  of  phase  and  hence  yielding  a  low  correla¬ 
tion  at  zero  time  lag.  In  contrast  to  a  scalar  analysis,  the  vector  analysis  combines  coherent  U  and 
V  winds  into  a  single  EOF  mode.  The  first  three  spatial  eigenvectors  (Figure  4)  are  rotated  so 
that  the  temporal  expansion  of  each  mode  is  aligned  along  its  principal  axis.  Mode  1  winds 
describe  63.4%  of  the  total  observed  low  frequency  variability,  are  aligned  along  the  gulf  axis, 
and  are  spatially  coherent  over  the  study  region  accounting  for  at  least  50%  of  the  variance  at 
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each  station.  These  winds  are  most  energetic  at  the  island  stations  W2  and  W5,  away  from  coa¬ 
stal  boundaries.  Mode  2  (19.5%)  isolates  the  strong  cross-gulf  winds  observed  at  W2.  Modes  3 
(11.3%)  and  4  (4.4%,  not  shown)  represent  winds  at  Puerto  Penasco  and  Guaymas  respectively 
that  are  probably  associated  with  localized  gales. 

In  later  comparisons  with  moored  observations,  the  dominant  component  of  mode  1 
wind  is  used  to  represent  the  along-gulf  wind  forcing.  A  representative  amplitude  of  this  forcing 
at  the  mooring  transect,  however,  is  not  obvious  (Figure  4)  given  that  winds  appear  to  be  twice  as 
strong  at  Isla  Tortuga  (W5)  than  at  the  coast  (W4,  W6).  This  also  implies  the  presence  of  a 
cross-gulf  wind  curl  which  could  have  a  significant  effect  on  gulf  circulation.  One  possible 
explanation  for  the  varying  wind  speeds  is  the  difference  in  station  elevations  at  the  transect.  W5 
was  at  the  top  of  Isla  Tortuga  at  approximately  150-m  elevation,  W4  was  on  a  100-m  bluff  at  the 
coast,  while  W6  was  on  a  30-m  high  island  less  than  1  km  from  the  coast.  It  is  therefore  unclear 
whether  PAMII  winds  are  representative  of  winds  near  the  water  surface.  To  investigate  these 
issues,  PAMII  winds  are  compared  with  winds  measured  30  m  above  the  water  surface  during  14 
overflights  of  the  transect  in  March  1984.  Overflight  winds  are  averaged  into  1-km  sections 
along  the  cross-gulf  track  and  EOFs  obtained  from  the  spatially-averaged  covariance  matrix. 
This  analysis  complements  the  wind  station  EOFs;  whereas  the  PAMII  data  consists  of  3  spatial 
measurements  at  the  transect  for  over  300  days,  the  overflight  data  has  140  spatial  measurements 
on  12  different  days.  The  overflight  and  PAMII  mode  1  spatial  patterns  (Figure  5)  are  normalized 
to  unit  value  at  Isla  Tortuga  (W5),  leaving  the  temporal  functions  in  physical  units.  The 
overflight  pattern  shows  that  winds  are  much  more  uniform  across  the  gulf  than  implied  by 
PAMII  winds.  The  mode  1  temporal  behavior  of  PAMII  and  overflight  along-gulf  winds  are  in 
good  agreement  indicating  that,  despite  the  poor  temporal  resolution  of  the  overflights  and  the 
possibility  of  diurnal  wind  aliasing,  the  same  variability  is  isolated  in  each  EOF.  In  addition,  the 
comparable  amplitudes  imply  that  PAMII  winds  at  Isla  Tortuga  are  similar  to  winds  near  the  sur¬ 
face.  The  overflight  winds  indicate  then  that  along-gulf  wind  amplitude  can  be  considered  uni¬ 
form  across  the  transect  with  negligible  cross-gulf  shears,  and  that  this  amplitude  is  best 
represented  by  the  Isla  Tortuga  PAMII  station  rather  than  the  coastal  stations. 

An  important  property  of  the  along-gulf  wind  is  its  relationship  to  the  surface  pressure 
gradient.  Geostrophic  winds,  computed  from  the  difference  in  atmospheric  pressure  between 
Guaymas  and  Santa  Rosalia,  are  in  excellent  agreement  (correlation  =  0.87)  with  PAMII  mode  1 
winds  (Figure  6).  Geostrophic  winds  are  stronger  than  observed  winds  with  a  regression  of  0.7  at 
Isla  Tortuga,  where  mode  1  winds  are  at  a  maximum.  This  is  consistent  with  observations  made 
over  water  in  the  German  Bight  (Hasse  and  Wagner,  1971),  where  the  regression  between  surface 
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and  geostrophic  winds  was  0.6.  The  geostrophic  wind  can  be  simply  obtained  from  pre-existing 
weather  stations  near  Guaymas  and  Santa  Rosalia  and  is  used  in  later  comparisons  with  moored 
observations  in  the  absence  of  measured  winds. 

Diurnal  winds  dominate  higher  frequency  wind  fluctuations  with  peak  speeds  ranging 
from  3  to  5  m/s.  These  winds  are  oriented  perpendicular  to  shore  with  winds  measured  at  the 
coast  10-20%  stronger  than  at  the  island  stations.  C.  E.  Dorman  (personal  communication) 
describes  the  diurnal  circulation  as  a  conventional  sea  breeze  along  the  mainland  coast  and  an 
upslope  valley  effect  along  the  Baja  California  coast  driven  by  the  heating  of  mountain  slopes. 
No  significant  difference,  however,  in  the  magnitude  of  diurnal  winds  at  opposite  sides  of  the  gulf 
is  observed. 

Temperature 

Near-surface  temperatures  in  the  Gulf  of  California  exhibit  a  large  seasonal  range  with 
especially  high  summer  temperatures  in  comparison  to  Pacific  waters  at  the  same  latitude  across 
the  Baja  California  peninsula  (Robinson,  1973).  The  temperature  range  measured  10  m  below 
the  surface  is  about  18°C  at  Guaymas  and  13°C  at  Santa  Rosalia  (Figure  7).  The  1983  El  Nino 
event  (Robles  and  Marinone.  1987)  is  evident  as  Guaymas  temperature  is  approximately  5°C 
warmer  in  the  winter-spring  of  1983  than  in  the  same  period  of  the  following  year.  Although  the 
averaged  near-surface  temperature  over  the  period  May  1983  to  May  1984  is  about  the  same 
(23°C)  at  Guaymas  and  Santa  Rosalia,  the  range  is  larger  at  Guaymas  with  warmer  summer  and 
cooler  winter  temperatures.  This  same  cross-gulf  difference  is  also  found  in  the  island  region 
north  of  the  transect  P7  -  P10  (Figure  1).  The  larger  temperature  range  on  the  mainland  coast  is 
consistent  with  seasonal  upwelling  (Thorade,  1909).  Upgulf  summer  winds  are  up/downwelling 
favorable  along  the  Baja  Califomia/mainland  coast  of  the  gulf  resulting  in  the  cooler  tempera¬ 
tures  at  Santa  Rosalia.  The  situation  is  reversed  in  winter  when  downgulf  winds  result  in  cooler 
temperatures  on  the  mainland  coast.  The  tendency  then  is  for  seasonal  winds  to  moderate  the 
seasonal  solar  heating  of  gulf  water  on  the  Baja  California  shelf  while  enhancing  the  extremes  on 
the  mainland  shelf. 

The  occurrence  of  seasonal  up/downwelling  is  also  apparent  in  the  contrasting  vertical 
temperature  profiles  on  opposite  sides  of  the  gulf  (Figure  8).  During  the  early  summer,  when 
thermistor  chains  were  operating  on  both  sides,  temperatures  below  100  m  are  slightly  warmer  at 
Santa  Rosalia  (1°C)  while  above  this  depth  Guaymas  temperatures  are  much  warmer  (5°C)  with  a 
well-developed  thcrmocline  region  between  40-70-m  depth  (0.2°C/m).  In  constrast,  the  strongest 
gradient  at  Santa  Rosalia  (0.17°C/m)  occurs  in  the  upper  20  m.  During  the  winter  the  vertical 


structure  has  reversed  across  the  gulf.  At  Guaymas,  average  temperatures  for  November  1984 
resemble  the  summer  upwelling  profile  at  Santa  Rosalia  with  the  absence  of  a  seasonal  thermo- 
cline  region.  While  the  thermistor  chain  at  Santa  Rosalia  was  not  operational  during  the  winter,  a 
CTD  cast  from  November  1984  (PC6003  from  Bray  et  al.,  1986),  near  the  shelf  break,  has  a 
deeper  thermocline  region  consistent  with  downwelling  favorable  winds  for  the  Baja  California 
shelf  during  this  time  of  year. 

Temperature  fluctuations  about  the  seasonal  cycle  are  considerable  at  Guaymas  during 
the  summer,  particularly  at  low  frequencies  (Figure  9).  In  contrast,  temperature  fluctuations  at 
Santa  Rosalia  during  the  early  summer  are  not  as  large.  On  both  shelves,  the  largest  fluctuations 
(Figure  10)  occur  at  the  depths  of  the  highest  mean  vertical  temperature  gradient  which,  as  noted 
above,  is  stronger  and  deeper  in  the  water  column  at  Guaymas  in  summer.  When  scaled  by  the 
mean  temperature  gradient  at  each  depth,  low  frequency  fluctuations  are  more  uniform  with  depth 
on  both  shelves.  Assuming  that  temperature  fluctuations  are  primarily  associated  with  the  verti¬ 
cal  displacement  of  isotherms,  the  weighted  profile  provides  an  approximate  measure  of  typical 
isotherm  displacements.  At  Santa  Rosalia,  these  displacements  increase  from  5  m  near  the  sur¬ 
face  to  10  m  below  80-m  depth  while  values  at  Guaymas  range  from  10-15  m. 

The  spatial  and  temporal  behavior  of  Guaymas  low  frequency  temperature  fluctuations, 
minus  the  seasonal  cycle,  are  described  with  EOFs  for  the  main  1984  deployment.  Mode  1  (65% 
of  the  variance)  describes  temperature  fluctuations  which  occur  in  phase  over  the  entire  shelf  and 
are  strongest  in  the  thermocline  region  (Figure  11).  There  is  a  significant  correlation  (0.75'> 
between  the  temporal  behavior  of  mode  1  and  Guaymas  adjusted  sea  level,  such  that  a  rise  in  coa¬ 
stal  sea  level  is  associated  with  a  sinking  of  isotherms  over  the  shelf  (Figure  1 1).  A  regression  of 
sea  level  with  the  depth  of  the  21°  isotherm  at  the  Guaymas  shelf  break,  which  is  generally  below 
the  surface  mixed  layer  for  most  of  the  summer,  shows  that  a  10-cm  rise  in  sea  level  corresponds 
to  a  drop  of  approximately  20  m  for  this  isotherm.  The  same  is  true  for  isotherms  deeper  in  the 
water  column.  In  addition,  this  rek, ions  hip  between  sea  level  and  temperature  occurs  in  winter 
as  well  when  temperature  70  m  below  the  surface  at  the  Guaymas  midshelf  (M7)  and  adjusted  sea 
level  are  also  correlated  (0.62).  Mode  2  (15.5%)  has  a  spatial  function  with  one  zero  crossing 
with  depth  such  that  fluctuations  approximately  above  and  below  the  thermocline  (-30  m)  are 
180°  out  of  phase  (Figure  1 1).  The  temporal  expansion  for  this  mode  is  significantly  correlated 
(0.65)  with  bottom  pressure  at  the  100-m  isobath  (Figure  11).  Visually,  the  similarity  between 
the  two  time  series  is  strongest  during  a  pressure  rise  in  July  which  will  be  seen  in  the  following 
section  to  be  the  largest  propagating  event  of  the  1984  summer.  This  large  event  is  also  evident 
in  mode  1  indicating  that  temperature  fluctuations  associated  with  propagating  events  are  not 
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isolaud  by  a  single  mode.  The  signal  of  the  propagating  events  below  40-m  depth  in  the  first  two 
modes  sum  constructively  indicating  that  sub-thermocline  isotherm  displacements  are 
emphasized  during  these  events. 

A  similar  EOF  analysis  is  performed  on  residual  low  frequency  temperature  on  the  Santa 
Rosalia  shelf  (Figure  12)  revealing  a  less  coherent  spatial  structure  than  at  Guaymas.  Like  Guay- 
mas,  low  frequency  temperature  changes  appear  to  be  uniform  across  the  shelf  as  fluctuations  at 
the  100-m  (Mil)  and  200-m  (M10)  moorings  are  similar  at  depth.  Mode  1,  however,  describes 
only  53%  of  the  total  observed  variability  and  has  a  zero  crossing  with  the  depth  in  contrast  to 
mode  1  at  Guaymas.  The  temporal  behavior  of  mode  1  temperature  is  only  weakly  correlated 
(0.39)  with  adjusted  sea  level  at  Santa  Rosalia  (Figure  12). 

On  both  shelves  as  well  as  in  the  basin,  low  frequency  winds  and  temperature  functions 
are  not  significantly  correlated  in  any  season.  Instead,  much  of  the  temperature  signal  below  the 
upper  mixed  layer  is  associated  with  propagating  variability.  This  is  somewhat  surprising  given 
the  evidence  for  seasonal  upwelling  presented  above  and  the  reported  upwelling  events  seen  in 
satellite  images  (Badan-Dangon  et  al.,  1985).  One  possibility  is  that  the  seasonal  upwelling  sig¬ 
nal  may  be  a  result  of  upwelling  primarily  occurring  about  specific  topographic  features.  Upwel¬ 
ling  plumes  along  the  mainland  coast  often  occur  north  of  the  Guaymas  moorings  (Badan- 
Dangon  et  al.,  1985),  closer  to  the  complicated  topography  of  the  island  region.  Upwelled  water 
may  then  be  mixed  and  advected  to  the  transect  array  on  longer  time  scales,  creating  the  observed 
seasonal  signal  without  a  correlation  between  wind  and  temperature  at  low  frequencies.  The  pro¬ 
pagating  signal  in  temperature  is  apparent  in  maximum  correlations  (Table  2),  with  Topolobampo 
leading  Guaymas  by  40  hours  and  Guaymas  leading  Santa  Rosalia  by  8.5  days.  The  nature  of  this 
long  lag  across  the  gulf  is  discussed  below  in  the  context  of  coastal-trapped  waves  (section  6). 

Standard  deviations  of  temperature  fluctuations  in  the  diumal/inertial  (0.92-1.08  cpd) 
and  semi-diurnal  (1.92-2.08  cpd)  bands  for  Guaymas  and  Santa  Rosalia  during  the  summer  1984 
deployment  (Figure  13)  show  that,  as  for  lower  frequencies,  the  largest  fluctuations  occur  at 
depths  of  strongest  stratification  on  both  shelves.  When  scaled  by  the  mean  temperature  differ¬ 
ence  at  each  depth,  displacements  range  from  2-8  m  on  each  shelf  with  weak  diumal/inertial 
fluctuadons  at  Guaymas  (Figure  13).  Semi-diurnal  band  EOFs  at  Guaymas  reveal  that  much 
(83%)  of  the  temperature  signal  is  nearly  in  phase  at  all  depths  with  a  vertical  structure  similar  to 
the  rms  profile  shown  in  Figure  13.  The  temporal  behavior  of  this  mode  has  the  same  spring- 
neap  modulation  as  the  semi-diurnal  band  surface  tide  with  sea  level  leading  temperature  by 
approximately  8  hours.  These  fluctuations  are  larger  by  a  factor  of  5  than  the  signal  which  could 
be  produced  by  the  barotropic  dde  moving  the  thermistor  chain  vertically  by  the  mean 
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temperature  gradient.  Instead  these  fluctuations  are  likely  due  to  internal  tides.  The  observations 
described  here  are  not  adequate  to  determine  the  origin  of  this  signal.  Semi-diurnal  fluctuations 
on  the  Santa  Rosalia  shelf,  however,  are  incoherent  with  the  surface  tide  and  do  not  exhibit  the 
same  coherent  vertical  structure  as  at  Guaymas.  The  same  is  true  of  tidal  band  fluctuations  in  the 
basin  and  on  the  mainland  shelf  at  Topolobampo.  It  is  unclear  why  this  signal  is  only  apparent  at 
Guaymas.  Diurnal  fluctuations  on  both  shelves  are  more  stochastic  and  appear  unrelated  to  the 
barotropic  tide,  diurnal  winds,  or  daily  heating  and  cooling. 

Bottom  Pressure  and  Sea  Level 

Sea  level  has  a  seasonal  range  of  approximately  0.5  m  in  the  gulf,  reaching  a  maximum 
in  July  and  a  minimum  in  February.  The  same  range  is  not  observed  in  bottom  pressure  below 
100  m  suggesting  that  seasonal  changes  in  sea  level  are  associated  with  the  heating  of  upper  layer 
waters.  To  verify  this,  changes  in  sea  surface  height  about  a  mean  value  is  estimated  from  the 
200-m  thermistor  chain  and  bottom  pressure  sensor  (P5)  at  the  Guaymas  shelf  break.  Density 
changes  in  the  gulf  are  roughly  correlated  with  changes  in  temperature  (Bray,  1988).  CTD 
profiles  near  the  Guaymas  shelf  provide  a  means  of  estimating  density  directly  from  temperature, 

p(z)  =  1.028  -  0.854  x  10-4  T(z)  -  0.355  x  10'6  7(z)2  (1) 

_c 

with  a  standard  error  of  2.3  x  10  g/cm  .  Pressure  near  the  surface  is  then  computed  from  the 
hydrostatic  equation  using  bottom  pressure  and  density  from  the  thermistor  chain, 

P (z0)  =  P (zb)  - g  j  p(z)dz  (2) 


where  zb  =  200 m ,  and  related  to  sea  level. 


(3) 


Over  the  duration  of  the  thermistor  chain  deployment,  the  range  in  T)  is  approximately  0.2  m  in 
both  observed  and  computed  sea  level.  In  equations  2  and  3,  P  ( zb )  and  Pmm  contribute  only  a 
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few  cm  to  the  seasonal  range  in  T)  confirming  that  density  changes  in  the  upper  200  m  account  for 
most  of  the  observed  seasonal  sea  level  signal. 

At  subtidal  frequencies,  the  poleward  propagating  events  described  by  Christensen  et  al. 
(1983)  and  Enfield  and  Allen  (1983)  are  evident  in  sea  level  and  bottom  pressure  records  along 
the  mainland  coast  during  this  experiment  (Figure  14).  Atmospheric  pressure  fluctuations 
account  for  no  more  than  10%  of  sub-tidal  sea  level  variability  along  the  coast  and  cannot  explain 
the  propagating  signal.  Local  winds  are  also  uncorrelated  with  sea  level  and  pressure  in  the  gulf. 
Most  of  the  events  are  present  as  far  south  as  Salina  Cruz  (about  2000  km  south  of  Topolobampo) 
and  travel  approximately  non-dispersively  into  the  Gulf  of  California.  Sea  level  at  Santa  Rosalia 
does  not  show  the  event  propagating  down  the  Baja  coast.  Travel  times  and  hurricane  tracks  sug¬ 
gest  that  events  at  Cabo  San  Lucas  (CB)  are  not  waves  continuing  down  the  Baja  California  pen¬ 
insula  but  are  instead  forced  by  the  hurricanes  themselves  as  they  skirt  the  southern  tip  of  Baja 
California.  Phase  speeds  are  estimated  from  neighboring  sea  level  and  pressure  sensors  using 
lagged  correlations  and  times  between  peak  events  (Table  3)  with  the  events  shown  in  Figure  14. 
The  estimates  arc  in  good  agreement  with  previous  values  obtained  by  Enfield  and  Allen,  1983. 

On  the  mainland  shelf,  sub-tidal  bottom  pressure  fluctuations  decay  rapidly  offshore. 
Bottom  pressure  fluctuations  at  the  Guaymas  200-m  shelf  break  (P5),  approximately  12  km  from 
shore,  are  weak  compared  to  coastal  adjusted  sea  level  particularly  during  the  large  propagating 
events  of  July  1984  (Figure  15).  In  contrast,  pressure  fluctuations  at  100-m  depth  (P6)  and 
adjusted  sea  level  at  Santa  Rosalia  are  similar  in  magnitude.  The  Guaymas  P5  bottom  pressure, 
however,  does  not  take  into  account  the  contribution  of  density  on  pressure,  or  internal  pressure 
changes.  Therefore,  pressure  just  below  the  sea  surface  and  at  100-m  depth  above  the  bottom 
pressure  sensor  P5  is  computed  from  equation  2  (zQ  =  0  m  and  100  m)  using  the  200-m  thermis¬ 
tor  chain  at  the  Guaymas  shelf.  When  compared  at  similar  depths,  pressure  fluctuations  are 
nearly  identical  across  the  shelf  (Figure  15).  The  contribution  of  sea  level  changes  to  pressure  is 
compensated  by  adjustments  in  the  density  field  such  that  pressure  fluctuations  at  100  m  are 
reduced  to  half  the  near  surface  amplitudes,  and  200-m  pressure  is  nearly  constant.  That  the  spa¬ 
tial  scale  of  variation  of  bottom  pressure  is  not  the  same  as  for  sea  surface  elevation  is  evident  in 
the  first  EOF  (90%  of  total  variance)  of  the  the  sub-tidal  pressure  field  at  Guaymas  (Figure  16). 
The  horizontal  e-folding  scale  is  on  the  order  of  50  km,  compared  to  the  30-km  deformation 
radius  given  by  Christensen  et  al.  (1983),  while  the  same  decay  occurs  in  160  m  in  the  vertical. 
On  the  Santa  Rosalia  shelf,  density  changes  do  not  significantly  contribute  to  fluctuations  in  the 
pressure  field  in  the  upper  100  m. 
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Although  local  winds  are  uncorrelated  with  sea  level  and  pressure  at  all  locations  around 
the  gulf,  the  wind  response  may  be  masked  by  the  more  energetic  propagating  variability.  In 
order  to  account  for  this,  near  surface  pressure  at  Isla  Tiburon  (P10)  is  estimated  from  adjusted 
sea  level  at  Guaymas  and  along-gulf  wind  stress  in  a  multiple  regression.  The  Guaymas  time 
series  is  lagged  in  time  to  obtain  maximum  correlation  with  sea  level  at  Isla  Tiburon  and  to 
approximate  the  portion  of  the  signal  associated  with  travelling  waves.  Wind  stress  is  obtained 
from  mode  1  winds  described  above  using  the  relation  given  by  Large  and  Pond  (1981).  Ideally 
if  there  is  a  residual  relationship  between  wind  and  sea  level,  the  multiple  regression  should 
return  a  significant  coefficient  for  winds.  The  regression  explains  81%  of  the  variance  at  Isla 
Tiburon,  but  nearly  90%  of  the  explained  variance  is  attributed  to  lagged  Guaymas  adjusted  sea 
level.  The  regression  coefficient  for  along-gulf  wind  stress  is  small  and  not  significantly  different 
than  zero.  This  result  of  no  local  wind  forcing  is  similar  for  all  pressure  records  around  the  gulf 
and  for  different  lags  of  winds  relative  to  pressure. 

Currents 

Representative  subinertial  currents  on  the  Guaymas  and  Santa  Rosalia  shelves  are  com¬ 
pared  in  Figure  17  for  the  entire  two-year  experiment  Currents  at  Santa  Rosalia  were  measured 
10  m  below  the  surface  at  the  90-m  isobath  while  the  Guaymas  time  series  is  a  compilation  of 
currents  measured  closest  to  the  100-m  isobath  (maximum  separation  =  3  km).  The  alongshelf 
directions  are  chosen  as  the  principal  axes  for  each  time  series  which  agree  well  with  the  orienta¬ 
tion  of  local  isobaths  on  each  shelf.  Currents  are  strongly  polarized  with  80-90%  of  the  total  vari¬ 
ability  directed  in  the  alongshelf  direction.  Alongshelf  currents  are  both  stronger  and  more  vari¬ 
able  at  Guaymas  than  at  Santa  Rosalia  with  peak  speeds  of  75  and  50  cm/s,  and  rms  speeds  of  15 
and  10  cm/s  repectively.  Current  spectra  at  Guaymas,  Santa  Rosalia,  and  the  Guaymas  Basin,  are 
red  with  significant  peaks  at  the  tidal  harmonics  (Figure  18).  Energy  levels  are  higher  in  the  sum¬ 
mer  than  winter  over  all  frequencies,  although  a  direct  comparison  is  only  valid  at  Santa  Rosalia 
given  the  different  locations  of  Guaymas  observations.  Subinertial  currents  are  most  energetic  at 
Guaymas  in  both  seasons  while  higher  frequency  currents  are  more  similar  in  energy  content  at 
the  three  locations.  The  difference  in  high  frequency  energy  between  Guaymas  and  Santa  Rosalia 
during  the  winter  is  probably  due  to  stronger  currents  in  1984  relative  to  1983  rather  than  a  differ¬ 
ence  between  the  two  shelves. 

Alongshelf  currents  at  Santa  Rosalia  (Figure  17)  appear  to  have  a  slight  seasonal  signal 
with  a  tendency  for  downgulf  flows  during  summer  and  upgulf  during  winter.  This  is  opposite  to 
the  direction  of  seasonal  winds.  To  investigate  seasonal  patterns  in  the  current  field,  mean 
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currents  are  computed  separately  for  the  summer  (June-  September),  when  winds  are  blowing  to 
the  north,  and  the  remainder  of  the  year,  when  winds  are  blowing  southward.  Currents  during 
these  two  periods  are  considered  to  have  a  significant  seasonal  mean  if  the  mean  differs  from  zero 

by  at  least  3  standard  errors,  with  the  standard  error  of  a  sample  mean  defined  as  (<v  2  >  /  N  )0  5 

/ 

where  v  is  the  standard  deviation  of  the  flow  and  N  is  the  number  of  independent  samples 
(Davis,  1976).  Along-gulf  currents  have  significant  mean  flows  in  the  basin  at  50-m  depth  (order 
10  cm/s)  and  on  the  Santa  Rosalia  shelf  at  10-m  and  30-m  depth  (5-8  cm/s)  in  both  summer  and 
winter,  with  50-m  basin  currents  flowing  in  the  direction  of  seasonal  winds  and  currents  at  Santa 
Rosalia  directed  against  the  wind.  This  occurrence  of  cyclonic  summer  and  anti-cyclonic  winter 
surface  circulation  in  the  Guaymas  Basins  was  observed  by  Bray  (1988)  using  hydrographic  data. 
This  pattern  does  not  appear  to  extend  onto  the  Guaymas  shelf  where  mean  flows  are  not 
significant.  A  causal  relationship  between  seasonal  winds  and  the  sense  of  rotation  of  Guaymas 
Basin  circulation  is  not  known. 

At  low  frequencies,  currents  vary  considerably  across  the  Guaymas-Santa  Rosalia  tran¬ 
sect  (Figure  19).  On  the  Guaymas  shelf,  surface  currents  decay  rapidly  in  amplitude  towards  the 
coast  with  rms  speeds  decreasing  from  27  cm/s  at  the  shelf  break  (M8)  to  10  cm/s  at  midshelf 
(M6)  and  2  cm/s  at  the  inner-shelf  (M4).  Similarly  the  alongshelf  polarization  of  currents 
increases  towards  the  coast  with  the  percentage  of  alongshelf  to  total  rms  current  speed  increasing 
from  74%  at  the  shelf  break  to  86%  near  the  coast.  Correlations  of  alongshelf  currents  between 
the  three  moorings  are  significant  but  low  (Table  4).  Maximum  correlations  occur  with 
alongshelf  currents  at  the  midshelf  (M6)  leading  currents  at  the  shelf  break  (M8)  by  32  hours  and 
inner  shelf  (M4)  by  52  hours.  The  currents  at  these  locations  are  most  similar  during  several 
energetic  events  in  late  June  and  July,  in  part  associated  with  propagating  events  described  below. 
While  the  vertical  structure  of  the  current  field  is  not  well  sampled,  the  alongshelf  currents  at  M4 
and  in  the  upper  30  m  at  M6  are  correlated  at  different  depths  (Table  4)  and  are  similar  in  ampli¬ 
tude.  In  addition,  alongshelf  currents  at  10  m  and  70  m  in  100-m  total  depth  (M7),  measured  dur¬ 
ing  winter  and  late  summer  of  1983,  arc  also  well  correlated  with  comparable  amplitudes  indicat¬ 
ing  that  alongshelf  flow  is  reasonably  barotropic  over  the  shelf.  In  contrast,  correlations  of 
cross-shelf  flows  are  typically  not  significant  between  the  various  sensors. 

Current  measurements  in  the  basin  (M9)  are  most  energetic  near  the  surface  with  rms 
speeds  of  22  cm/s  at  50-m  depth  compared  to  7  cm/s  at  300  m  and  5  cm/s  at  850  m.  Currents  at 
50  m  are  directed  primarily  along  the  gulf  and  are  not  correlated  with  deeper  flows.  Currents  at 
300  m  and  500  m  are  nearly  identical  with  principal  axes  directed  more  in  the  across  gulf  direc¬ 
tion,  primarily  due  to  one  strong  cross-gulf  flow  event  in  July.  Weak  currents  at  850  m  appear 


14 


unrelated  to  flows  higher  in  the  water  column. 

Currents  on  the  steep  Santa  Rosalia  shelf  are  more  spatially  uniform  than  at  Guaymas 
(Figure  19)  with  rms  speeds  ranging  from  10-15  cm/s.  Currents  seperated  by  2  km  across  the 
shelf  (M10  and  Ml  1)  and  10  m  below  the  surface  are  correlated  (Table  4)  in  both  the  alongshelf 
(0.94)  and  cross-shelf  (0.71)  directions.  The  vertical  structure  of  currents  at  Ml  1  appears  barotro- 
pic  in  Figure  19  particularly  during  the  energetic  flow  in  late  July.  The  correlation  is  larger 
between  flow  at  10  m  and  30  m,  than  flows  between  10  m  and  70  m,  due  in  part  to  the  change  to 
upgulf  surface  flow  in  late  fall  which  does  not  occur  at  70-m  depth.  The  weaker  alongshelf  flow 
observed  in  the  summer/fall  of  1983  (5-9  cm/s  rms  speeds)  is  even  less  uniform  with  depth  as 
correlations  between  flow  at  10  m  and  65  m  were  not  significant. 

Santa  Rosalia  and  Guaymas  low  frequency  currents  are  uncorrelated  in  all  seasons 
(Table  4),  confirming  the  visual  impression  (Figures  17  and  19)  that  shelf  flows  on  opposite  sides 
of  the  gulf  are  unrelated.  On  the  mainland  coast,  alongshelf  currents  at  Topolobampo  and  Guay¬ 
mas  are  correlated  during  both  summer  deployments  with  Topolobampo  leading  by  about  29 
hours,  similar  to  the  lags  found  in  temperature  and  pressure  correlations.  Winter  currents  at 
Topolobampo  were  not  available  for  comparison.  Although  currents  are  unrelated  across  the 
gulf,  basin  flow  is  correlated  with  alongshelf  currents  at  both  Santa  Rosalia  and  Guaymas. 
Alongshelf  currents  in  the  basin  at  50-m  depth  (M9  Figure  20)  have  longer  time  scales  (4.8  days) 
than  at  10-m  depth  on  the  Guaymas  shelf  (3.1  days).  The  significant  negative  correlation  (-0.54) 
is  most  apparent  in  June  and  July  during  4  periods  of  strong  opposite  flowing  currents.  Wind  and 
current  fluctuations  appear  related  at  this  time  although  a  causal  relationship  is  unclear  given  that 
wind  events  are  not  particularly  energetic  and  later  periods  of  strong  downgulf  winds  do  not  have 
an  associated  current  response.  These  ,  sporadic  occurrences  of  opposing  currents  may  indicate 
eddy-like  flow  features  near  the  Guaymas  shelf  break.  Although  satellite  images  during  these 
particular  events  are  unavailable,  the  presence  of  eddies  in  this  area  is  supported  by  satellite 
images  showing  entrainment  patterns  in  sea  surface  temperature  with  scales  on  the  order  of  30  km 
(C.  Paden,  personal  communication).  The  correlation  between  alongshelf  currents  at  Santa 
Rosalia  (10-m  depth  at  Ml  1)  and  at  the  basin  is  also  negative  (-0.61).  The  correlation  on  this  side 
of  the  gulf  is  due  primarily  to  opposing  flows  in  July  and  early  fall.  The  opposing  flow  is  remark¬ 
ably  steady  over  the  month  of  October  and  is  related  to  the  seasonal  pattern  described  above. 

Currents  above  the  inertial  frequency  have  typical  rms  speeds  of  2-3  cm/s  in  this  region 
of  the  gulf.  Current  spectra  (Figure  18)  show  significant  peaks  at  the  diurnal  and  semi-diurnal 
bands  at  Santa  Rosalia,  while  at  Guaymas  only  the  semi-diurnal  band  is  signficant.  Like  low  fre¬ 
quency  currents,  tidal  flows  are  polarized  in  the  alongshelf  direction  with  more  cross-shelf  energy 


15 


further  from  the  coast.  Tidal  band  currents  are  not  coherent  at  different  sensor  locations  on  either 
shelf  or  in  the  basin.  Surprisingly,  this  is  true  at  Guaymas  where  semi-diurnal  temperature 
fluctuations  suggest  the  presence  of  Internal  tides.  Alongshore  current  at  the  Guaymas  shelfbreak 
(10-m  depth  M8)  is  coherent  (0.74)  with  the  mode  1  semi-diurnal  band  temperature  signal  with 
currents  lagging  by  roughly  2  hours. 


5.  MOMENTUM  TERMS 
Cross -She  If 

A  number  of  recent  field  observations  have  verified  that  a  near  geostrophic  balance 
exists  at  the  coast  between  alongshelf  current  and  cross-snelf  pressure  gradient, 

v  -(f  pTldP  Idx  (4) 

The  first  verifications  using  direct  measurements  of  bottom  pressure  were  made  during  the  Nan¬ 
tucket  Shoals  Flux  Experiment  (NSFE,  Brown  et  al.,  1985)  and  the  Coastal  Ocean  Dynamics 
Experiments  (CODE,  Brown  et  al.,  1987,  Winant  et  al.,  1987).  The  CODE  experiments  were 
made  along  a  straight  coast  in  an  upwelling  regime  characterized  by  strong  local  wind  forcing, 
weak  stratification,  and  nearly  barotropic  flows.  NSFE  was  conducted  on  a  wide  shelf,  10  times 
the  width  of  the  Guaymas  shelf,  with  moderate  local  wind  forcing  and  high  summer  stratifications 
similar  to  those  found  in  the  gulf.  In  both  of  these  studies,  suitably  averaged  alongshelf  current 
observations  were  in  excellent  agreement  with  the  geostrophic  current  derived  from  the  cross¬ 
shelf  gradient  in  bottom  pressure  measurement..  The  depth-averaged,  or  barotropic  transport, 
was  compared  to  geostrophic  current  for  the  CODE  observations,  while  cross-shelf  averages  at 
various  depths  were  compared  for  the  NSFE  observations. 

The  geostrophic  current  at  Guaymas  is  obtained  from  the  pressure  difference  between 
adjusted  coastal  sea  level  and  surface  pressure  at  the  shelfbreak  (P5,  see  Figure  15).  The  pres¬ 
sure  gradient  is  obtained  in  a  similar  manner  on  the  Santa  Rosalia  shelf  with  surface  pressure 
computed  above  the  100-m  pressure  sensor  (P6)  using  temperature  observations  from  the  three 
current  meters.  This  correction  for  internal  pressure  changes,  however,  is  negligible  at  Santa 
Rosalia.  In  this  construction  of  the  pressure  gradient,  the  geostrophic  current  is  a  measure  of  the 
transport  near  the  surface.  Observed  current  at  Guaymas  is  taken  as  an  average  of  the  alongshelf 
currents  at  M8,  M6  and  M4,  10  m  below  the  surface.  Santa  Rosalia  current  is  from  the  10-m 
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sensor  at  Ml  1  in  100  m  total  water  depth.  Since  pressure  gradients  are  found  to  have  a  significant 
seasonal  component  that  is  associated  with  errors  in  computing  density  directly  from  temperature 
rather  than  with  a  physical  signal,  all  time  series  involved  in  pressure  gradient  comparisons  are 
band-passed  (0.03-0.33  cpd).  If  the  geostrophic  balance  is  valid,  the  observed  and  geostrophic 
current  at  the  sea  surface  should  be  well  correlated  with  a  regression  near  1.0.  Correlations 
between  observed  and  geostrophic  currents  (Figure  21)  on  both  sides  of  the  gulf  are  signficant 
(0.47  at  Guaymas,  0.41  at  Santa  Rosalia)  with  regression  coefficients  of  0.71  and  1.09  at  Guay- 
mas  and  Santa  Rosalia  respectively.  Although  the  time  series  are  only  moderately  correlated, 
many  of  the  flow  features  at  Guaymas,  particularly  during  the  large  propagating  event  of  July,  are 
well  represented  by  the  geostrophic  current.  A  comparison  of  the  geostrophic  current  with  the 
individual  current  records  at  Guaymas  reveals  that  the  highest  correlation,  0.44,  occurs  at  M8, 
decreasing  to  0.37  at  M6,  with  an  insignifcant  correlation  at  M4  indicating  that  geostrophic 
currents  are  not  uniform  across  the  shelf.  An  outer  shelf  geostrophic  current  maximum  was  also 
found  in  the  NSFE  observations.  The  correlation  at  Santa  Rosalia  does  not  change  significantly  if 
the  depth-averaged  flow  is  used  in  place  of  the  single  current  record.  When  the  geostrophic 
current  is  obtained  from  the  pressure  gradients  measured  along  the  bottom  as  in  the  aforemen¬ 
tioned  studies,  correlations  with  observed  alongshelf  currents  are  no  longer  significant  emphasiz¬ 
ing  the  need  to  account  for  internal  pressure  effects  in  the  present  data  set. 

In  a  further  test  of  geostrophy,  the  difference  in  adjusted  sea  level  across  the  gulf 
between  Guaymas  and  Santa  Rosalia  is  examined.  This  cross-gulf  pressure  gradient  is  compared 
with  the  unweighted  average  of  alongshelf  currents  at  Santa  Rosalia  (Mil),  the  basin  mooring 
(M9),  and  Guaymas  (M8).  Including  weaker  nearshore  flow  at  Guaymas  does  not  significantly 
alter  this  average.  The  geostrophic  current  from  the  cross-gulf  pressure  gradient  and  average 
current  from  the  three  current  records  are  well  correlated  (0.70)  indicating  that  the  net  surface 
flow  through  a  cross-gulf  transect  can  be  represented  by  a  simple  pressure  gradient  (Figure  22). 
This  result  gives  insight  into  several  aspects  of  gulf  flow.  First,  most  of  the  pressure  gradient  sig¬ 
nal  is  due  to  pressure  fluctuations  at  Guaymas  with  only  weak  fluctuations  at  Santa  Rosalia  (Fig¬ 
ure  15a).  The  pressure  signal  at  Guaymas  has  been  shown  to  consist  primarily  of  propagating 
variability  along  the  mainland  coast  with  very  little  local  wind  contribution.  Hence  most  net 
subinertial  transport  at  the  surface  past  the  Guaymas-Santa  Rosalia  transect  occurs  during  pro¬ 
pagating  sea  level  events.  Secondly,  most  of  the  net  transport  takes  place  at  or  near  the  Guaymas 
shelf.  The  portion  of  flow  on  the  shelves  correlated  with  the  basin  flow  (see  Figure  20)  cancels 
out  in  the  averaging  with  the  remaining  net  flow  comprised  mainly  of  Guaymas  current.  The 
correlation  of  the  cross-gulf  pressure  gradient  with  Guaymas  current  alone  is  (0.58).  Much  of  the 
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current  in  the  basin  and  at  Santa  Rosalia  must  be  then  associated  with  recirculating  flow  patterns 
not  contributing  to  net  surface  transport.  If  the  main  surface  transport  is  indeed  associated  with 
wave  events  at  Guaymas,  it  may  be  inferred  from  the  mode  1  characteristics  of  the  density  field  at 
Guaymas  that  compensating  flow  for  the  net  transport  pictured  in  Figure  22  takes  place  deeper  in 
the  water  column,  as  with  an  internal  Kelvin  wave.  Deep  currents  near  the  Guaymas  shelf  were 
not  available  for  confirmation.  Fourth,  the  correlation  is  remarkably  high  given  the  limited  sam¬ 
pling  of  the  current  field.  This  lends  confidence  to  the  assumption  that  the  major  patterns  of 
along-gulf  flow  have  been  sampled,  or  alternatively,  that  any  energetic  flows  not  adequately  sam¬ 
pled  are  associated  with  recirculating  cells. 


Alongshelf 

The  linear  alongshelf  momentum  balance  for  subinertial  flow  is: 

dv  ,  -1  dp  ^  Bxy 

—  +fu  =  — -j-  +  -3— 
dt  p0  dy  dz 


(5) 


where  the  principal  sources  of  ly  are  surface  wind  stress  and  bottom  friction.  The  relative  mag¬ 
nitudes  of  terms  in  this  equation  vary  over  different  shelf  regions  (Brink,  1987),  however  when 
local  wind  stress  is  strong,  winds  and  alongshelf  flow  are  generally  related.  Lentz  and  Winant 
(1986)  found  that  the  alongshelf  momentum  balance  for  depth-averaged  subinertial  flow  varies 
across  the  southern  California  shelf.  On  the  inner  shelf,  wind  stress  is  balanced  by  bottom  fric¬ 
tion  while  in  deeper  water,  in  the  absence  of  strong  local  winds,  alongshelf  currents  are  balanced 
by  the  alongshelf  pressure  gradient.  A  similar  analytic  result  was  obtained  by  Mitchum  and 
Clarke  (1986)  who  found  that  the  inner  shelf  balance  between  surface  wind  stress  and  bottom 
friction  extends  offshore  to  total  water  depths  of  approximately  2.5  times  the  Ekman  e-folding 
depth.  A  rigorous  examination  of  the  momentum  balance  in  the  gulf  is  not  possible  given  the 
sparse  mooring  array.  Instead,  the  relative  importance  of  local  winds  and  pressure  gradients  in 
forcing  alongshelf  flows  is  examined  on  both  sides  of  the  gulf. 

A  multivariate  regression  of  the  form 

dv  /dt  =  ^  dP  /dy  +  (32t>  (6) 


is  performed  on  each  current  record  at  the  transect  during  the  main  deployment.  The  alongshelf 
pressure  gradient  on  the  mainland  shelf  is  obtained  from  the  difference  between  adjusted  sea  level 
at  Guaymas  and  near  surface  pressure  at  Isla  Tiburon  (P10),  and  on  the  Baja  California  shelf 
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between  adjusted  sea  level  at  Santa  Rosalia  and  near-surface  pressure  at  San  Francisquito  (P7).  It 
is  noted  that  point  measurements  of  the  flow  are  compared  with  pressure  gradients  over  approxi¬ 
mately  300  km  distances.  Bottom  friction  is  not  considered  since  current  measurements  are  typi¬ 
cally  near  the  surface.  In  addition,  fu  is  not  included  given  the  noisy  records  and  poor  spatial 
coherence  of  cross-shelf  currents  over  either  shelf.  In  each  regression,  the  two  inputs  are  shifted 
in  time  to  obtain  maximum  correlation  with  dv  /dt .  The  hindcast  skill  equivalent  to  the  percen¬ 
tage  of  dv  /dt  variance  explained  in  each  regression,  is  significant  on  both  shelves  (Table  5), 
although  less  than  10%  of  the  variance  is  accounted  for  at  Santa  Rosalia  (Ml  1).  The  regression  is 
most  successful  in  explaining  currents  at  the  Guaymas  shelf  break  (M8)  accounting  for  over  40% 
of  the  variance.  On  the  Guaymas  shelf,  nearly  all  of  the  explained  variance  is  due  to  the 
alongshelf  pressure  gradient  with  winds  contributing  a  negligible  percentage.  The  negative 
values  of  are  consistent  with  equation  5.  The  pressure  gradient  lags  dv/dt  at  Guaymas  as 
would  be  expected  for  a  poleward  propagating  pressure  signal  as  the  pressure  difference,  taken 
between  Guaymas  and  Isla  Tiburon,  is  centered  north  of  Guaymas.  In  contrast,  the  pressure  gra¬ 
dient  along  the  Baja  California  peninsula  is  not  related  to  currents  at  Santa  Rosalia.  This  does  not 
necessarily  imply  that  the  pressure  gradient  in  unimportant  in  driving  flow  along  this  shelf.  The 
poor  relationship  may  indicate  that  pressure  gradients  over  shorter  length  scales  than  the  300  km 
sensor  seperation  are  important  between  Santa  Rosalia  and  San  Francisquito.  Winds  are  also 
unrelated  to  dv  /dt  on  the  Santa  Rosalia  shelf  and  in  the  basin.  The  negligible  effect  of  local 
wind  is  similar  if  v  is  used  instead  of  dv  /dt  in  equation  6.  This  is  true  on  the  Guaymas  inner 
shelf  where  a  relationship  between  wind  stress  and  bottom  friction  (assumed  proportional  to  v ) 
was  observed  in  the  aforementioned  studies.  The  lack  of  a  wind  response  also  occurs  whether 
local  winds  (W4  or  W6)  are  used  instead  of  the  spatially  uniform  EOF  wind,  or  if  different  lags 
between  wind  and  dv/dt  are  used  in  the  regression. 

The  reason  for  the  weak  correlation  between  winds  and  currents  is  examined  by  compar¬ 
ing  wind  magnitudes,  stratification,  and  current  wind  response  in  the  gulf  with  other  shelf  studies 
(Table  6).  Currents  in  the  first  two  studies  in  Table  6  were  characterized  as  having  a  strong  wind 
response,  as  indicated  by  typical  correlations  between  alongshelf  winds  and  currents  near  the  sur¬ 
face  at  a  midshelf  location  (Winant  et  al„  1987;  Beardsley  et  al.,  1985).  Common  to  both  studies 
was  a  combination  of  weak  stratification  and  strong  winds.  A  moderate  wind  response  was 
reported  in  the  next  two  studies.  On  the  southern  California  shelf,  winter  winds  were  weak  but 
the  lack  of  stratification  apparently  allows  a  measureable  wind  response  in  the  current  record 
(Lentz  and  Winant,  1986).  Conditions  at  NSFE  were  strongly  stratified  during  the  summer  but 
winds  were  also  fairly  strong  (Beardsley  et  al.,  1985).  The  remaining  studies  are  typified  by  a 
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weak  wind  response  and  in  most  cases  characterized  by  high  stratifications  and  weak  winds.  Dur¬ 
ing  the  CUEA  experiment  on  the  Peruvian  shelf  (Brink  et  al.,  1983),  stratification  was  similar  to 
CODE-2  but  weak,  fluctuating  winds  and  the  presence  of  coastal-trapped  wave  variability  led  to  a 
low  wind  response.  On  the  southern  California  shelf  (Lentz  and  Winant,  1987),  the  relationship 
between  winds  and  currents  was  inconclusive  although  a  typical  wind  response  was  not  observed, 
negative  correlations  occurred  during  summer.  In  the  context  of  these  studies,  the  absence  of  a 
measurable  wind-driven  subinertial  circulation  in  the  Gulf  of  California  is  not  unusual  given  the 
weak  winds  and  strong  stratification  found  throughout  most  of  the  year. 


6.  A  COASTAL-TRAPPED  WAVE  EVENT 

The  importance  of  remotely  forced,  propagating  variability  along  the  mainland  shelf  of 
the  Gulf  of  California  is  evident  in  both  lagged  correlations  of  temperature,  pressure,  and  current 
observations  separated  along  the  coast  (section  4),  and  in  the  dynamic  relationship  between  the 
alongshore  pressure  gradient  and  current  accelerations  at  Guaymas  (section  5).  In  this  section, 
moored  observations  during  the  largest  wave  of  the  summer  of  1984  are  described,  with  the  goal 
of  establishing  the  qualitative  structure  and  dynamic  importance  of  these  events.  A  more  detailed 
comparison  of  observations  with  coastal-trapped  wave  theory  will  be  presented  in  a  later  paper. 

The  event  is  actually  a  sequence  of  two  20-cm  sea  level  elevations  (Figure  14)  in  July 
1984  associated  with  hurricanes  Genevieve  and  Fausto  which  were  active  south  of  the  gulf  region 
(Gunther,  1985).  As  noted  by  Christensen  et  al.  (1983),  events  often  consist  of  2-3  peaks  in  sea 
level,  attributable  to  the  grouped  occurrence  of  storms  (Gray,  1982).  The  poleward  propagation 
of  this  event  is  evident  in  sea  level  and  bottom  pressure  records  along  the  mainland  coast  as  far 
south  as  Salina  Cruz.  The  event  propagates  non-dispersively  from  Topolobampo  to  Isla  Tiburon 
with  an  amplitude  of  20  cm  and  an  alongshore  scale  on  the  order  of  the  length  of  the  gulf.  The 
apparent  damped  amplitude  at  Topolobampo  is  due  to  the  depth  of  the  sensor  (i.e.,  no  adjustment 
for  internal  pressure,  Figure  15).  The  event  appears  to  continue  into  the  northern  gulf  as  is  evi¬ 
dent  in  the  record  at  San  Francisquito  although  now  somewhat  altered  in  temporal  shape.  Sea 
level  at  Santa  Rosalia  does  not  show  the  event  propagating  down  the  Baja  California  coast.  A 
rise  in  sea  level  of  5  cm  occurs  at  Santa  Rosalia  approximately  1  to  2  days  after  the  leading  edge 
of  the  event  reaches  Guaymas  and  lasts  for  the  duration  of  the  event.  This  occurs  during  an  event 
of  similar  amplitude  in  1983,  and  was  also  observed  by  Christensen  et  al.  (1983,  see  their  Figure 
8)  for  a  1976  wave  event.  This  may  imply  that  larger  events  do  span  the  width  of  the  gulf. 
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The  temperature  and  flow  behavior  at  the  Guaymas-Santa  Rosalia  transect  during  this 
wave  event  is  illustrated  in  Figure  24.  Alongshelf  currents  at  Topolobampo  (Ml,  100-m  isobath) 
are  directed  up  the  gulf  with  amplitudes  of  roughly  50  cm/s,  and  are  barotropic.  At  Guaymas,  the 
wave  event  occurs  after  a  period  of  strongly  fluctuating  current,  particularly  at  the  shelf  break 
(see  M8  in  Figure  20).  Nevertheless  currents  at  the  mid  (M6)  and  outer  (M8)  shelf  have  the  same 
upgulf  flow.  There  is  a  slight  phase  shift  with  the  alongshore  current  at  M6  leading  M8  by 
roughly  10  hours.  The  strongest  nearshore  (M4)  currents  of  the  summer  occur  during  the  event 
but  are  directed  downgulf  counter  to  offshore  currents  and  the  geostrophic  pressure  gradient.  One 
explanation  may  be  some  type  of  recirculation  cell  associated  with  the  rough  bathymetry  at  the 
coast.  The  M4  mooring  was  located  just  southward  of  a  small  cape  and  offshore  island.  Deeper 
shelf  currents  were  not  available  although  the  barotropic  alongshelf  current  structure  observed  at 
Topolobampo  was  also  observed  at  Guaymas  (M7  at  the  100m  isobath)  during  a  simlar  event  in 
1983.  In  the  basin  (M9),  there  is  a  current  away  from  the  mainland  shelf  at  300  m  and  500  m 
when  the  event  passes  Guaymas  although  the  double-peaked  nature  of  the  event  seen  on  the 
Guaymas  shelf  is  not  evident  in  the  basin.  No  changes  in  alongshelf  current  or  temperature  wave 
response  is  observed  in  the  basin. 

At  Santa  Rosalia  (Mil),  the  most  energetic  currents  of  the  summer  (50  cm/s)  occur 
approximately  36  hours  after  the  large  event  has  passed  the  array  at  Guaymas.  This  current, 
directed  downgulf,  does  not  have  the  same  temporal  shape  as  the  wave  event  nor  is  it  apparent  in 
the  alongshelf  pressure  gradient  between  Santa  Rosalia  and  San  Francisquito.  The  same  down¬ 
gulf  current  response  occurs,  however,  during  an  event  of  similar  amplitude  in  the  summer  of 
1983.  The  fact  that  the  current  is  not  associated  with  a  pressure  gradient  between  the  island 
region  and  Santa  Rosalia  suggests  that  the  current  is  not  directly  related  to  either  the  wave  event 
or  with  a  buildup  of  a  sea  level  head  as  the  wave  enters  the  shallow  northern  gulf.  One  possibility 
is  a  cyclonic  circulation  cell  setup  in  the  basin  by  strong  upgulf  currents  on  the  mainland  shelf 
during  the  passage  of  energetic  wave  events  (Niels  Christensen,  personal  communication).  Such 
a  structure  was  observed  spanning  the  width  of  the  gulf  near  Topolobampo  during  a  drogue  study 
in  August  1978  (Emilsson  and  Alatorre,  1980).  An  energetic  event  was  present  in  the  gulf 
approximately  two  weeks  before  the  study  was  made.  Downgulf  currents  at  Santa  Rosalia  after 
the  July  1984  event  persist  for  approximately  two  weeks.  The  existence  of  such  a  closed  flow 
pattern  would  constitute  one  of  the  major  exchange  mechanisms  of  waters  across  the  gulf  during 
the  summer  and  warrants  further  investigation. 

The  temperature,  or  equivalently  density,  behavior  observed  on  the  mainland  shelf  is 
similar  to  a  first  mode  internal  Kelvin  wave;  isotherms  across  the  shelf  sink  by  approximately  60 
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m  during  the  20-cm  rise  in  sea  level.  The  displacements  are  not  limited  to  the  thermocline  region 
but  occur  with  similar  amplitudes  to  at  least  130  m.  Temperature  does  not  change  noticeably  in 
the  upper  mixed  layer  at  10-m  depth.  Alongshelf  currents  at  Guaymas  (M61Q  and  M810  m)  are 
correlated  (-0.50)  with  mode  1  temperature  at  Guaymas,  pictured  in  Figure  11a,  over  the  entire 
deployment.  Tiiis  is  the  only  signficant  temperature-current  correlation  at  the  Guaymas-Santa 
Rosalia  transect  observed  during  this  study.  The  negative  correlation  is  consistent  with  the  first 
mode  internal  Kelvin  wave  as  a  poleward  alongshelf  current  corresponds  to  a  deepening  of  isoth¬ 
erms  over  the  shelf. 

Pressure  at  Guaymas,  as  well  as  temperature  (Table  3),  is  correlated  with  near  bottom 
temperatures  at  Santa  Rosalia  with  Guaymas  leading  Santa  Rosalia  temperature  by  approximately 
8  days  in  both  the  1983  and  1984  summer  deployments  (Figure  25).  In  1983,  the  energetic  events 
in  August  and  September  at  Guaymas  are  clearly  related  to  temperature  fluctuations  at  Santa 
Rosalia.  In  1984,  the  comparison  is  less  clear  although  the  large  July  event  is  apparent  on  both 
sides  of  the  gulf.  Temperature  fluctuations  at  all  depths  at  Santa  Rosalia  were  stronger  during  the 
summer  of  1984  than  1983  suggesting  that  the  propagating  signal  is  present  but  masked  by  other 
variability.  A  preliminary  study  of  coastal-trapped  wave  phase  speeds  in  the  gulf,  using  the  algo¬ 
rithm  of  Brink  and  Chapman,  1987,  shows  that  mode  1  waves  have  a  travel  time  of  roughly  8 
days  between  Guaymas  and  Santa  Rosalia  assuming  the  events  travel  around  the  head  of  the  gulf 
without  correction  for  coastline  curvature  or  the  blocking  effect  of  the  islands.  If  strong  pro¬ 
pagating  events  do  in  fact  travel  around  the  perimeter  of  the  gulf  to  Santa  Rosalia,  their  charac¬ 
teristics  appear  substantially  altered  with  only  weak  sub-thermocline  fluctuations  apparent  in  our 
observations  at  this  long  lag  time. 


7.  SUMMARY 

Moored  observations  indicate  that  local  winds  are  relatively  unimportant  as  a  direct 
forcing  mechanism  for  subinertial  shelf  circulation  in  the  Gulf  of  California.  Along-gulf  winds, 
which  are  to  a  large  extent  uniform  over  the  northern  gulf,  are  not  correlated  with  fluctuations  in 
the  current,  temperature,  or  pressure  fields  at  low  frequencies.  The  weak  and  variable  winds  and 
high  stratification  in  the  gulf  are  typical  of  regions  of  poor  wind  response. 

On  seasonal  time  scales,  differences  in  temperature  and  vertical  temperature  profiles 
across  the  gulf  reflect  the  occurrence  of  seasonal  upwelling.  Up/do wnwelling  occurs  in  the  sum¬ 
mer  on  the  Santa  Rosalia/Guaymas  shelf  with  the  opposite  situation  occurring  in  winter. 
Scasonal-avcragcd  mean  currents  are  weak  (<  10  cm/s).  Near  the  surface,  mean  currents  arc  in 
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the  same  direction  as  seasonal  wind  in  the  basin  as  reported  by  Bray  (1988),  and  against  the  wind 
on  the  Santa  Rosalia  shelf.  Mean  currents  at  Guaymas  are  not  significant. 

Low-frequency  current,  temperature,  and  pressure  fluctuations  are  stronger  on  the  main¬ 
land  than  Baja  California  shelf  and  are  not  generally  correlated  across  the  gulf.  Maximum  lagged 
correlations  on  the  mainland  coast  indicate  poleward  propagation  of  alongshelf  current,  tempera¬ 
ture,  and  pressure  signals.  Current  speeds  are  strongest  at  the  shelf  break  at  Guaymas  with  very 
weak  nearshore  flows.  The  pressure  field  at  Guaymas  displays  an  isostatic  relationship  between 
adjusted  sea  level  and  temperature  fluctuations  such  that  bottom  pressure  at  200-m  depth  is  nearly 
constant.  Adjusted  sea  level  and  temperature  fluctuations  in  the  upper  water  column  are  weakly 
correlated  on  the  Santa  Rosalia  shelf.  Alongshelf  currents  at  Guaymas  and  Santa  Rosalia  are 
correlated  with  the  cross-shelf  pressure  gradient  with  the  correlation  at  Guaymas  strongest  for 
currents  near  the  shelf  break.  Baroclinic  pressure  must  be  taken  into  account  when  determining 
geostrophic  currents  from  pressure  at  different  depths.  In  the  alongshelf  momentum  balance,  the 
time  derivative  of  alongshelf  currents  and  pressure  gradients  along  the  coast  are  correlated  at 
Guaymas  indicating  that  a  signifcant  percentage  of  the  current  variabilty  along  the  mainland  coast 
is  associated  with  nonlocal  forcing.  Correlations  between  currents  and  pressure  gradients  are 
again  highest  for  currents  at  the  shelf  break.  This  relationship  between  currents  and  alongshelf 
pressure  gradients  was  not  observed  at  Santa  Rosalia  where  propagating  variability  is  not 
observed  in  sea  level  and  bottom  pressure  records. 

Tidal  and  inertial  currents  are  weak  (order  5  cm/s)  in  this  region  of  the  gulf.  Semi¬ 
diurnal  temperature  fluctuations  on  the  Guaymas  shelf  have  a  zero  mode  vertical  structure  and  are 
coherent  with  sea  level  displaying  in  amplitude  with  the  spring-neap  cycle.  Alongshelf  currents 
are  weakly  coherent  with  this  signal.  These  fluctuations  are  believed  to  represent  internal  tides 
although  they  are  only  observed  on  the  Guaymas  shelf. 

Storm-generated  coastal  waves  have  phase  speeds  on  the  order  of  200-300  km/day  as 
reported  by  Enfield  and  Allen  (1983).  Along  the  mainland  shelf  during  the  passage  of  a  20-cm 
event,  the  alongshelf  current  on  the  mid  shelf  is  roughly  50  cm/s  and  is  barotropic  in  structure. 
Currents  near  the  surface  are  strongest  at  the  shelf  break  and  decay  rapidly  towards  the  coast  with 
slight  phase  differences  across  the  shelf.  The  density  field  responds  like  a  first  mode  internal 
baroclinic  wave  with  a  20-cm  rise  in  sea  level  associated  with  a  40-m  depression  of  isopycnals 
throughout  much  of  the  water  column.  The  same  current  and  temperature  signals  are  not  found  at 
Santa  Rosalia  during  coastal-trapped  wave  events.  Alongshelf  currents  are  energetic  after  large 
events  but  are  not  associated  with  observed  alongshelf  pressure  gradients.  Sub-thermocline  tem¬ 
peratures  at  Santa  Rosalia  and  Guaymas  are  correlated  with  Guaymas  leading  Santa  Rosalia  by 
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7-9  days  which  is  consistent  with  the  computed  travel  time  of  these  waves  around  the  head  of  the 
gulf. 

Our  conclusion  is  that  shelf  circulation  on  either  side  of  the  gulf  is  largely  independent 
of  each  other.  Basin  flow  is  related  to  flow  patterns  on  both  the  Guaymas  and  Santa  Rosalia 
shelves.  Guaymas  and  the  basin  share  energetic  counter  flows  near  the  shelf  break  with  order  1 
week  time  scales,  while  longer  period  flows  at  seasonal  time  scales  run  counter  at  Santa  Rosalia 
and  the  basin.  The  mainland  shelf  appears  to  derive  much  of  its  energy  from  disturbances  which 
propagate  in  from  the  Pacific  ocean.  This  is  certainly  true  during  large  storm  generated  events, 
but  also  appears  to  be  the  case  for  weaker  variability  judging  from  the  good  comparison  between 
currents  and  alongshelf  pressure  gradients.  During  large  storm-generated  events,  a  response  is 
detected  on  the  Baja  shelf  although  it  is  different  in  character  than  on  the  mainland  shelf.  The 
cross-gulf  net  geostrophic  current  is  correlated  with  a  simple  average  of  along-gulf  currents 
between  Guaymas  and  Santa  Rosalia.  Most  of  this  net  transport  occurs  on  the  Guaymas  shelf 
during  propagating  events. 

It  is  unclear  to  what  extent  shelf  circulation  on  the  mainland  coast  of  the  Gulf  of  Califor¬ 
nia  is  influenced  by  remote  wind  forcing  other  than  summer  storm  events.  Shelf  measurements 
along  the  eastern  Pacific  coast  at  Washington  (Battisti  and  Hickey,  1984),  Oregon  (Cutchin  and 
Smith,  1973) ,  California  (Chapman,  1987),  and  Peru  (Smith,  1978)  all  document  the  importance 
of  remotely  forced  variability.  Enfield  and  Allen  (1983)  also  detected  slower  poleward  propaga¬ 
tion  over  shorter  length  scales  in  winter  when  tropical  storms  were  not  active.  Shelf  circulation 
in  the  Gulf  of  California  may  be  unique  from  other  marginal  seas  of  similar  proportions  in  that  it 
has  a  rather  long  coastal-trapped  wave  fetch  along  the  Mexican  Pacific  coast  from  which  to  draw 
energy.  One  might  imagine  a  very  different  current  climate  if  for  example  the  gulf  were  open  to 
the  north  and  closed  to  direct  poleward  propagation.  The  role  of  winds  south  of  the  central  gulf 
needs  to  be  considered  in  future  studies  of  shelf  circulation  in  this  region. 
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TABLE  1.  Deployment  record  for  instrumer’ation  from  1982-1984. 


Site 

Water  Depth 

(m) 

Instrument 

Depth/Elevation(+) 

(m) 

Deployment 

Length 

(days) 

Deployment  1 : 

1 1/82-5/83 

M7 

100 

10,  701 

164 

Deployment  2: 

Ml 

100 

102,  701 

153 

5/83-11/83 

M5 

35 

102 

130 

M7 

100 

102, 701, 75 

186 

Mil 

86 

10, 65 1 

190 

PI 

100 

100 

169 

P2 

35 

35 

187 

P6 

86 

86 

188 

Deployment  3: 

M7 

100 

102, 70* 

190 

1 1/83  5/84 

Mil 

90 

10, 651 

190 

W1 

- 

+20 

443 

W2 

- 

+50 

350 

W3 

- 

+0 

125 

W4 

- 

+  150 

424 

W5 

- 

+100 

446 

W6 

- 

+  10 

322 

Deployment  4: 

Ml 

100 

10,70 

107 

5/84-11/84 

M4 

15 

5.  12 

191 

M6 

50 

10,  20,  302 

188 

M8 

200 

10 

188 

M8  t -chain 

200 

10,  20,  30,  40,  50.  60, 70,  85.  100, 

115,  130,  145,160,175 

188 

M9 

1000 

50,  1002,  3001, 

5001, 8501 

179 

M10 

200 

10 

120 

MlOt-chain 

200 

10, 20, 30, 40. 50. 60, 70, 85, 100 

115, 130, 145, 160, 175 

58 

Mil 

90 

10,  30, 70 

188 

PI 

100 

100 

187 

P4 

100 

100 

381 

P5 

200 

200 

188 

P6 

90 

90 

187 

P7 

7 

7 

194 

P10 

7 

7 

191 

1  -  inclinometer 

2  -  temperature  only 
+  -  above  sea  level 


TABLE  2.  Correlations  for  low  frequency  temperature  records  during  the 
summer/fall  deployments  of  1983  and  1984.  Instrument  location  is 
represented  by  mooring  number  followed  by  depth  of  sensor  beneath  the 
surface.  Maximum  correlations,  if  different  than  zero-lag  correlation,  and 
time  lag  for  maximum  correlation  are  also  listed.  A  positive  lag  means 
that  temperature  at  the  first  sensor  leads  the  second. 


Instrument 

Maximum 

Lag 

Location 

Correlation 

Correlation 

(days) 

1983 


Topdobampo-Guaymas 


M1to-M7to 

0.78* 

0.0 

Ml70-M775 

M1,oo-M775 

0.52* 

0.62* 

1.38 

0.54* 

0.73* 

1.63 

Topolobampo-Santa  Rosalia 

M1to-MU.o 

Ml* ’MU* 

M1too-MU« 

0.66* 

0.0 

0.04 

0.67* 

8.42 

-0.02 

0.69* 

9.29 

Guaymas-Sanu  Rosalia 

M710-MI110 

0.50* 

0.56* 

7.13 

MVMII«6 

0.14 

0.62* 

7.50 

0.08 

0.71* 

7.33 

1984 


Topol  obampo-Guaymas 


M1to-M8to 

0.16 

0.0 

M17o-M870 

0.45* 

0.64* 

1.79 

^too-^too 

M1,oo-M8.oo 

0.47* 

0.67* 

1.58 

0.45* 

0.64* 

1.79 

Topolobampo-Santa  Rosalia 

MVMU70 

0.07 

0.39* 

9.42 

0.03 

0.46* 

9.21 

Guaymas-Basin 

M8„-M9„ 

-0.16 

0.0 

S 

s' 

1  f 
00 

s 

-0.21 

0.0 

Guaymas-Sanu  Rosalia 


M8io’MUio 

-0.08 

0.0 

MVMU30 

-0.21 

0.0 

M87o-MH7o 

-0.03 

0.55* 

8.96 

M7too-MI1*, 

0.03 

0.55* 

8.46 

Basin-Sanu  Rosalia 

M930  '  M11jo 

0.09 

0.0 

M9*  -MU* 

_ U2B _ 122 _ 

0.36 

0.0 

above  95%  significance  level 


TABLE  3.  Phase  speeds  of  propagating  pressure  signals  from  lagged  correlations  and  time  lags  between  peak 
events  for  May-September  1983  and  1984.  Values  obtained  by  Enfield  and  Allen  (1983)  for  the  1971  summer 
are  included  for  comparison. 
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1983 

1984 

1971 

Lag 

Correlation 

Peak 

Events 

Lag 

Correlation 

Peak 

Events 

Cross 

Spectra 

Peak 

Events 

AC-TO 

316 

344 

283 

313 

250 

316 

TO-GY 

204 

212 

252 

224 

225 

229 

GY-IT 

260 

252 

AC  =  Acapulco 
TO  =  Topolobampo 
GY  =  Guaymas 
IT  =  Is  la  Tibonon 
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TABLE  4.  Correlation!,  as  in  TABLE  2,  for  alongshelf  currents  during  the  4  instrument  deployments. 


Location 

Instrument 

Correlation 

Maximum 

Correlation 

Lag  (days) 

11/82-5/83 

GY 

M7  -  M7___ 

0.82* 

0.0 

lCkn  70m 

5/83-11/83 

GY-TO 

ML.  -  M7 

0.49* 

0.60* 

1.1 

GY-SR 

M7--M11- 

70m  65m 

-0.28 

0.0 

11/83-5/84 

GY-SR 

M7_.  -Mil.. 

70m  10m 

0.07 

0.0 

5/84-11/84 

GY 

M4.  -  M6„__ 

-0.14 

0.51* 

-2.2 

M8j”  -  Ml- 

0.14  • 

0.24* 

1.0 

M4  -  M9  — 

-0.26* 

0.0 

M6*“  -  M§— 

20m  10m 

0.35* 

0.52* 

1.3 

SR 

MU,,  -Mil,. 

0.86* 

0.0 

MU "  -Mlt 

0.34* 

0.0 

Mil— -MU— 

0.55* 

0.0 

M10—  -Mil— 

0.94* 

0.0 

lOn  10m 

TO 

MI  -Ml 

0.78* 

0.0 

10m  70m 

BA 

M9  -  M9-„ 

0.90* 

0.0 

M9”--M9^" 

0.41* 

0.0 

300m  850m 

TO-GY 

Ml..  ■  M8 

10m  10m 

-0.08 

0.30* 

1.9 

GY-BA 

M620m-M950m 

-0.15 

0.46* 

3.3 

IChn  50m 

-0.48* 

0.53* 

1.0 

BA-SR 

M9„  -MU.. 

-0.54* 

0.65* 

2.7 

M9—  -Mil— 

■0.55* 

2.1 

_ dftn  _ 

*  -  above  95%  significance 
TO  -  Topolobampo 
GY  -  Guaymas 
BA  -  Guaymas  Basin 
SR  -  Santa  Rosalia 


TABLE  5.  Multivariate  analysis  of  the  alongshelf  pressure  gradient  (dP/dy)  and  along -gulf  wind  stress  regressed  on  dV/dt  at  vari¬ 
ous  moorings  from  equation  6.  Hindcast  skill  is  the  amount  of  variance  explained  in  the  regression,  %  variance  represents  the  vari¬ 
ance  explained  by  each  input. 
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dP/dv 

T » 

(xlOO) 

Hindcast 

Skill 

Pi 

Lag 

(hours) 

Vo 

Variance 

Pi 

Lag 

(hours) 

% 

Variance 

dV/dt 

M45 

11.2* 

-0.05  ±0.06 

0 

11.3 

0.02  ±0.03 

0 

1.0 

M620 

14.9* 

-0.27  ±0.27 

-37 

15.5 

-0.04±0.15 

0 

0.1 

M810 

41.6* 

-0.82  ±0.61 

-15 

33.0 

0.51  ±0.33 

0 

3.7 

M9jo 

9.1 

0.11  ±0.09 

0 

4.4 

0.06  ±0.06 

0 

4.9 

M,1to 

7.7* 

-0.09  ±0.05 

32 

4.3 

0.06  ±0.06 

0 

3.5 

M1130 

2.7* 

-0.05  ±0.03 

33 

2.6 

0.01  ±0.05 

0 

0.1 

M1170 

6.5 

0.00  ±0.03 

33 

2.3 

-0.06  ±0.05 

0 

5.3 

above  95%  significance  level 


TABLE  6.  Comparison  of  stratification,  wind  stress,  and  wind-current  correlations  measured  during  different 
field  studies. 


Experiment 

Brunt-Vaisala 

Frequency 

N  (cpd) 

<r>> 

Mean 

Aiongshelf 

Wind 

(PA) 

T»' 

Sl-Dev.  of 
Aiongshelf 
Wind 
(PA) 

(V,T>) 

Correlation  Between 
Near-surface,  Midshelf 
Aiongshelf  Current 
and  Wind 

CODE-2 

100 

-0.12 

0.12 

0.7  -  0.8 

NSFE79 

(summer) 

30 

0.01 

0.27 

0.63 

So.  California 
(winter) 

50 

- 

0.01 

0.41-0.69 

NSFE79 

(summer) 

2S0 

0.01 

0.05 

0.36 

Gulf  of  Calif, 
(winter) 

140-170 

-0.01 

0.02 

0.2 

Peru 
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Figure  1. 


(a)  The  Gulf  of  California  showing  mooring  (M),  bottom  pressure  (P),  and  PAMII 
meteorological  station  (W)  locations;  (b)  Locations  of  sea  level  observations. 
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Figure  3.  Low  frequency  vector  time  series  of  PAMII  wind  observations  showing  the  seasonal  wind 
pattern  over  the  Gulf  of  California.  Vectors  pointing  up  the  page  represent  winds  blowing 
up  the  Gulf  axis  (330°). 
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LOW  FREQUENCY  PRINCIPAL  AXES 


NODE  1  <63.450 


Figure  4. 


Principal  axes  for  low  frequency  wind  fluctuations  and  the  three  lowest  spatial  EOF  modes 
of  the  low  frequency  vector  wind  observations  during  the  time  period  March  1984- 
February  1985.  Station  W3  is  not  included  in  the  EOF  analysis.  Vector  amplitudes  are  in 
arbitrary  units.  Values  next  to  arrows  represent  the  percentage  of  variance  explained  in 
each  mode  at  that  particular  station.  Vectors  are  rotated  by  a  constant  value  for  each  mode 
to  coincide  with  the  principal  axis  for  the  temporal  expansion  of  that  mode. 


MODE  1  TEMPORAL 


de  1  EOFs  obtained  from  the  PAMII  (solid  line)  and  overflight  (triangles)  observations 
Lhe  Guaymas-Santa  Rosalia  transect  during  March  1984.  The  spatial  amplitudes  have 
:n  normalized  to  unit  value  at  Isla  Tortuga  (W5).  Temporal  amplitudes  are  in  physical 
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Figure  6.  Along-gulf  wind  at  W5  (thin)  compared  with  the  geostrophic  wind  obtained  from  the  sur¬ 
face  atmospheric  pressure  gradient  between  W4  and  W6. 
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Figure  7.  Subinertial  temperature  measured  10  m  below  the  surface  at  mid-shelf  locations  on  the 
Guaymas  and  Santa  Rosalia  shelves  and  the  difference  between  the  two  series. 
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Figure  9. 

Hourly-averaged  temperature  versus  depth  at  the  200-m  thermistor  chains  at  Guaymas 
(M8)  and  Santa  Rosalia  (MIO).  Temperature  series  are  equally  spaced  in  the  vertical. 
Mean  temperature  at  each  depth  is  pictured  in  Figure  8.  Mean  temperatures  at  M8  do  not 
change  significantly  when  averaged  over  the  entire  deployment. 
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Figure  10. 

(a)  Standard  deviations  of  low  frequency  temperature  fluctuations  versus  depth  at  the  Santa 
Rosalia  (dash)  and  Guaymas  (solid)  200-m  thermistor  chains  for  the  time  periods  pictured 
in  Figure  8.  (b)  The  same  profiles  scaled  by  the  average  vertical  temperature  gradient 
computed  at  each  depth  . 
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Figure  11. 

(a)  The  two  lowest  spatial  EOF  modes  of  low  frequency  temperature  at  the  Guaymas  shelf 
during  the  summer/fall  1984  deployment.  Values  are  proportional  to  the  amount  of  vari¬ 
ance  explained.  An  annual  harmonic  was  removed  from  each  time  series  before  forming 
the  unweighted  covariance  matrix  for  this  analysis.  Asterisks  indicate  sensor  positions! 

(b)  A  comparison  of  the  mode  1  temporal  expansion  and  Guaymas  adjusted  sea  level,  and 
the  mode  2  temporal  expansion  and  100-m  bottom  pressure  (P4).  Units  of  temporal 
expansions  are  arbitrary. 


46 


1984 


11  MAY  8  JUN 


i  !  ■ 

6  JUL 


i  1  i  '  _ I _ i — : — : — I — - 

3  AUG  31  AUG  28  SEP 


j _ L_i 

26  OCT 


Figure  12. 

The  two  lowest  spatial  EOF  modes  of  low  frequency  temperature  at  the  Santa  Rosalia  shelf 
during  May-July  1984  for  the  thermistor  chain  at  M10  (solid)  and  the  VMCM  mooring  at 
Mil  (dash).  Asterisks  indicate  sensor  locations.  Values  are  proportional  to  the  amount  of 
variance  explained  at  each  sensor,  (b)  A  comparison  of  the  mode  1  temporal  expansion 
and  Santa  Rosalia  adjusted  sea  level,  and  the  mode  2  temporal  expansion.  Temporal 
expansions  are  in  arbitrary  units. 


Figure  13. 

Standard  deviations  of  diurnal  (asterisks)  and  semi-diunal  (triangle)  band  temperature 
fluctuations  versus  depth  at  the  Santa  Rosalia  (dash)  and  Guaymas  (solid)  200-m  thermis¬ 
tor  chains  for  the  time  periods  pictured  in  Figure  8.  The  diurnal  band  is  defined  as  0.92- 
1.08  cpd  and  the  semi-diurnal  as  1.92-2.08  cpd.  (b)  The  same  profiles  scaled  by  the  aver¬ 
age  vertical  temperature  gradient  at  each  depth. 
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gure  5.  Adjusted  sea  level  and  bottom  pressure  at  100  m  (P4)  and  200  m  (P5)  at  Guaymas, 
adjusted  sea  level  and  90-m  bottom  pressure  (P6)  at  Santa  Rosalia,  and  Guaymas  com¬ 
puted  pressure  at  0  m  and  100  m  above  the  200-m  (PS)  sensor,  using  equation  2 
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Figure  17.  Low  frequency  currents  at  Santa  Rosalia  (Mll]0  )  and  Guaymas  (^7^  m  from  11/82- 
4/83.  M7„  5/83-8/83.  M7  9/83-10/83.  m  1 1/83-4/84,  average  Between  M62Q  m 

and  M8in  V  5/84-1 1/84). 
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Figure  19. 

Low  frequency  vector  time  series  of  currents  across  the  Guaymas-Santa  Rosalia  transect 
during  the  summer-fall  1984  deployment  and  mode  1  wind  stress.  A  vector  pointing  up 
the  page  corresponds  to  upgulf  flow. 


Figure  20.  Comparison  of  alongshelf  currents  in  the  basin,  on  the  Guaymas  and  Santa  Rosalia 
shelves,  and  mode  1  wind  stress  during  the  summer-fall  1984  deployment. 
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Figure  24. 

Current  and  temperature  observations  in  the  Gulf  of  California  during  the  most  energetic 
propagating  event  of  1984  as  seen  in  adjusted  set 1  vel  at  Guaymas. 
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Figure  25. 

Comparison  of  shallow  bottom  pressure  and  adjusted  sea  level  at  Guaymas  with  tempera¬ 
ture  at  Santa  Rosalia  for  the  summers  of  1983  and  1984.  Events  which  appear  to 
correspond  in  both  records  are  connected. 


CHAPTER  2 


A  COMPARISON  OF  PRESSURE,  CURRENT,  AND  DENSITY  OBSERVATIONS 
WITH  COASTAL-TRAPPED  WAVE  THEORY 

1.  INTRODUCTION 

Sub-tidal  sea  level  events,  with  positive  elevations  of  10-30  cm  and  10  day  time  scales, 
have  been  reported  propagating  poleward  along  the  Mexican  coast  and  into  the  Gulf  of  California 
(Christensen  et  al.,  1983;  Enfield  and  Allen,  1983,  henceforth  CDG  and  EA).  Based  on  estimates 
of  phase  speed,  offshore  extent,  and  mass  transport,  CDG  conclude  that  these  sea  level  signals  are 
the  surface  expression  of  hybrid  coastal-trapped  waves,  more  similar  to  baroclinic  Kelvin  waves 
than  barotropic  shelf  waves.  EA  use  cross-spectra,  case  history  and  multiple  regression  analysis 
to  show  that  these  waves  are  associated  with  tropical  storms  and  hurricanes  that  occur  during  the 
summer  and  early  fall  off  the  Pacific  coast  of  Mexico.  They  find  good  agreement  between  the 
observed  propagation  and  free  coastal-trapped  wave  phase  speeds  north  of  20°N,  while  faster  pro¬ 
pagation  south  of  20°N  is  consistent  with  the  speeds  of  the  traveling  storms.  These  waves  are 
somehow  dissipated  in  the  Gulf  of  California,  as  they  are  not  detected  in  tide  gauge  observations 
on  the  Pacific  Ocean  side  of  the  Baja  California  peninsula  (CDG). 

The  internal  characteristics  of  these  waves  have  been  observed  with  a  moored  array  in 
the  Gulf  of  California  (Merrifield  and  Winant,  1988,  or  thesis  chapter  1,  henceforth  MW).  On  the 
mainland  coast  during  an  energetic  event,  currents  are  strongest  at  the  shelf  break  and  decrease 
rapidly  towards  the  coast.  Consistent  with  the  hypothesis  that  the  events  are  wavelike,  MW  find 
the  alongshelf  current  is  in  approximate  geostrophic  balance  with  the  cross-shelf  pressure  gra¬ 
dient,  and  the  acceleration  of  the  alongshelf  current  is  well  correlated  with  the  alongshelf  pressure 
gradient  and  uncorrelatcd  with  local  winds.  Over  the  mainland  shelf,  a  sea  level  rise  of  10  cm 
corresponds  to  a  20-m  downward  displacement  of  isotherms  below  the  mixed  layer.  On  the  Baja 
California  shelf,  subthermocline  temperature  fluctuations  arc  signficantly  correlated  with  sea 
level  at  Guaymas  at  an  8  to  9  day  lag  suggesting  that  some  form  of  wave  energy  propagates 
around  the  head  of  the  gulf  and  returns  down  the  peninsula.  Sea  level  and  current  observations, 
however,  do  not  indicate  cquatorward  wave  propagation  on  this  side  of  the  gulf. 

In  this  chapter,  we  attempt  to  determine  if  this  propagating  variability  is  consistent  with 
a  forced  coastal-trapped  wave  model  (Clarke  and  Van  Gorder,  1986).  Due  to  the  lack  of  reliable 
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wind  observations  along  this  portion  of  the  Mexican  coast,  a  complete  hindcast  comparison  is  not 
possible.  Instead,  the  nature  of  the  more  energetic  storm-related  events  described  above  arc 
investigated  by  extrapolating  storm  winds  to  the  coast.  An  empirically  derived  propagating  sig¬ 
nal,  obtained  from  observed  coastal  pressure  and  sea  level,  is  extracted  from  the  moored  observa¬ 
tions  and  compared  to  the  model.  These  gulf  observations  provide  a  test  of  the  coastal-trapped 
wave  model  under  several  unique  conditions:  summer  stratification  is  very  strong  (buoyancy  fre¬ 
quency  =  250  cpd),  wave  events  tend  to  be  isolated  occurrences  that  can  be  tracked  over  several 
thousand  kilometers,  and  there  is  little  evidence  of  local  wind  forcing  in  the  gulf.  The  paper  is 
organized  as  follows.  A  brief  description  of  the  observations  and  the  coastal  wind  estimate  dur¬ 
ing  storms  is  given  in  section  2,  a  discussion  of  the  coastal-trapped  wave  model  and  the  free  wave 
parameters  for  this  region  is  contained  in  section  3,  the  empirical  propagating  signal  is  derived 
and  compared  to  model  hindcasts  in  section  4,  and  a  discussion  and  summary  follow  in  section  5. 


2.  OBSERVATIONS 

Moored  observations  of  currents,  temperature,  and  bottom  pressure  have  been  collected 
in  the  Gulf  of  California  during  a  joint  field  study  by  Centro  de  Investigation  Cientifica  y  dc 
Education  Superior  de  Ensenada  (CICESE)  and  the  Scripps  Institution  of  Oceanography  (SIO)  in 
1983-84  (Fig.  1).  A  cross-gulf  transect  between  Guaymas  on  the  mainland  coast  and  Santa 
Rosalia  on  the  Baja  California  peninsula  is  the  primary  study  site.  Current  observations  are  made 
with  vector  measuring  current  meters  (VMCMs),  while  temperature  is  measured  at  each  VMCM 
and  pressure  sensor  and  with  thermistor  chains  in  200  m  of  water  on  the  Guaymas  and  Santa 
Rosalia  shelves.  Winds  and  atmospheric  pressure  arc  measured  with  Portable  Automated 
Mesonet  (PAMII)  weather  stations  positioned  at  six  island  and  coastal  sites.  In  addition,  sea  level 
is  obtained  from  tide  gauges  maintained  by  CICESE  and  the  National  Autonomous  University  of 
Mexico  (UNAM),  and  atmospheric  pressure  south  of  the  gulf  is  obtained  from  the  World 
Meteorological  Organization.  A  complete  description  of  the  gulf  observations  is  contained  in 
MW. 

Reliable  coastal  wind  measurements  in  the  region  south  of  the  Guaymas-Santa  Rosalia 
transect  are  not  available  during  this  experiment.  Winds  measured  at  airports  near  the  coast  con¬ 
tain  numerous  gaps  and  spurious  points  and  are  not  considered  to  be  representative  of  winds  over 
water.  Moreover,  storm  events  are  typically  missing  from  these  wind  observations.  Large  scale 
gcostrophic  winds,  obtained  from  the  Fleet  Meteorological  Center,  do  not  have  the  spatial  resolu¬ 
tion  needed  to  detect  hurricane  winds  and  are  of  dubious  value  at  the  low  latitudes  in  this  study. 
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Thus,  in  order  to  study  the  coastal-trapped  wave  model  response  during  energetic  tropical  storm 
events,  storm  track  and  maximum  wind  speed  data,  obtained  from  the  National  Climatic  Data 
Center,  are  used  to  estimate  winds  at  the  coast.  During  a  tropical  storm  or  hurricane,  an  estimate 
of  storm  position  and  wind  speed  is  made  every  six  hours  from  satellite  data  (Dvorak,  1973),  with 
occasional  ship  and  aircraft  verification.  Given  the  approximate  storm  center  and  maximum  wind 
speed,  we  extrapolate  tangential  winds  to  the  coast  using  typical  radial  wind  speed  profiles  for 
hurricanes  and  tropical  storms  (Fig.  2)  given  by  Riehl  (1954).  Different  radial  profiles  for  tropi¬ 
cal  storm  force  winds  are  used  before  and  after  a  hurricane  under  the  assumption  that  a  storm  in 
formation  has  a  larger  spatial  extent  than  a  decaying  storm.  In  this  case,  we  use  the  hurricane 
radius  of  maximum  wind  speeds  after  the  disturbance  has  been  downgraded  to  a  tropical  storm. 
If  the  storm  docs  not  reach  hurricane  intensity,  however,  the  tropical  storm  profile  is  used  con¬ 
tinuously.  Storm  centers  that  are  greater  than  500  km  from  the  coast  and  tropical  depressions  arc 
not  included  in  the  analysis.  Wind  stress  is  then  computed  using  the  expression  of  Large  and 
Pond  (1981).  A  synthetic  time  scries,  which  only  accounts  for  storm  related  alongshore  winds,  is 
then  used  to  force  the  coastal  wave  model. 


3.  THE  COASTAL-TRAPPED  WAVE  MODEL 


The  dynamical  model  used  in  the  present  study  is  based  on  the  forced  barotropic  shelf 
wave  formulation  of  Gill  and  Schumann  (1974),  extended  to  include  stratification  (Gill  and 
Clarke,  1974),  cross-shelf  topography  (Clarke,  1977),  and  bottom  friction  (Brink  and  Allen,  1978; 
Brink,  1982;  Clarke  and  Van  Gorder,  1986).  Following  Clarke  and  Van  Gorder  (1986),  the 
linear,  long  wave,  wind-driven  response  at  the  coast  can  be  expressed  as  a  vorticity-type  equation 
for  pressure  P , 
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(1) 


where  /  is  the  coriolis  parameter,  N  (z)  is  the  Brunt-Vaisala  frequency,  and  x  and  z  arc  the 
cross-shelf  (positive  in  the  offshore  direction)  and  vertical  coordinates.  The  boundary  conditions 
are  no  flow  through  the  bottom,  a  rigid  lid  at  the  surface,  zero  depth-averaged  cross-shelf  flow,  u, 
at  a  depth  equal  to  three  times  the  Ekman  depth  (see  Mitchum  and  Clarke,  1986),  and  ux  =  0  at 
an  offshore  distance  L  which  is  typically  equal  to  twice  the  shclf/slope  width.  Equation  1  is 
solved  at  a  given  alongshclf  location  by  expanding  P  in  terms  of  the  complete  set  of  free  wave 
modes  Fj  (x ,  z ), 
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P(x,y,z,t)  =  X  Fj(x,z)tyj(y,t) , 
/= i 


(2) 


with  each  satisfying  a  forced,  first  order  wave  equation, 

oo 

Cf%  +  <$>jy  +  x °ij <>>  =bjty(y,t) , 
i=l 


(3) 


where  is  the  alongshelf  wind  stress,  Cj  is  the  free  wave  phase  speed,  bj  is  the  wind  coupling 
coefficient,  and  atj  are  the  frictional  parameters.  Friction  has  the  effect  of  coupling  the  modal 
solutions  in  (3),  as  well  as  providing  dissipation.  The  competing  effects  of  generation  by 
alongshelf  winds  and  frictional  dissipation  are  thus  integrated  as  the  wave  propagates  poleward. 
In  practice,  P  is  obtained  by  dividing  the  coastline  into  approximately  uniform  sections  (Fig.  3), 
and  for  each  section,  specifying  a  representative  bathymetry  and  stratification  and  computing  the 
free  wave  parameters  Fj (x  ,z),  Cj,  a and  bj  from  a  program  listed  in  Brink  and  Chapman 
(1985).  A  coastal  wall  is  placed  at  a  depth  of  10  m  to  approximate  the  coastal  boundary  condi¬ 
tion  of  Mitchum  and  Clarke  (1986a).  is  then  computed  at  a  given  alongshelf  location  from  (3) 
using  a  program  provided  by  Dr.  Allan  Clarke.  The  solution  for  §j(y  ,t)  requires  the 
specification  of  <]>;  (0,  r)  at  the  start  of  the  grid.  The  (0,  t)  used  in  this  study  are  described  in 
the  following  section.  Alongshelf  current  and  density  are  then  obtained  directly  from  P  using, 
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This  model  has  had  the  most  success  in  hindcasting  coastal  sea  level  and  alongshelf  currents 
(Battisti  and  Hickey,  1984;  Mitchum  and  Clarke,  1986b;  Chapman,  1987;  Church  et  al„  1986), 
although  current  speeds  are  usually  underprcdictcd.  Density  and  cross-shclf  current  hindcasts, 
however,  arc  generally  poor  (e.g.,  Chapman).  Chapman  (1987)  has  tested  the  sensitivity  of  the 
model  in  hindcasting  observations  made  during  the  CODE  experiment  off  northern  California. 
He  finds  the  model  is  most  sensitive  to  the  choice  of  wind  forcing  with  measured  coastal  winds 
yielding  better  hindcasts  than  the  larger  scale  gcostrophic  winds. 
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To  compute  the  free  wave  parameters,  bathymetry  charts  (reference)  are  used  to  obtain 
representative  cross-shelf  depth  profiles  for  each  section.  The  Brink  and  Chapman  program 
requires  the  specification  of  a  flat  bottom  beyond  the  shelf  slope  for  numerical  stability.  A  con¬ 
stant  depth  of  4000  m  is  used  for  sections  1  through  5  while  3500  m  is  used  for  sections  6  and  7 
as  the  ocean  floor  rises  near  the  mouth  of  the  gulf.  In  sections  8  through  12,  the  basin  shaped 
bathymetry  profile  of  the  Gulf  of  California  is  used  instead  of  a  flat  bottom.  The  boundary  condi¬ 
tions  are  no  flow  through  the  offshore  solid  boundary,  and  ux  =  0  at  x  =  L ,  where  L  is  the  width 
of  the  gulf.  The  rugged  bathymetry  of  the  island  region  of  the  northern  gulf  is  simplified  to  a 
smooth,  shallow  basin  of  approximately  200-m  depth.  Test  runs  confirm  that  the  modes  used  in 
this  study  are  stable  to  small  changes  in  the  gulf  bathymetry  profile.  A  further  restriction  of  the 
model  is  a  straight  coastline.  This  is  a  reasonable  assumption  along  most  of  the  Mexican  coast 
except  in  the  northern  gulf  where  the  coastline  orientation  changes  by  180°.  In  this  study,  the 
Gulf  of  California  will  simply  be  "unfolded"  with  the  coastline  assumed  straight  from  the  main¬ 
land  to  the  Baja  California  shelves.  Summer  stratification  profiles  for  sections  1  through  7  are 
obtained  from  historic  EASTROPAC  and  CALCOFI  hydrographic  data.  The  stratification  for 
sections  8  and  9  is  an  average  profile  obtained  from  the  thermistor  chain  on  the  Guaymas  shelf 
during  the  1984  summer,  sections  10  and  11  is  an  average  of  hydrographic  observations  made 
during  this  experiment  (Bray,  1987),  and  section  12  is  obtained  from  the  thermistor  and  mooring 
temperature  observations  on  the  Santa  Rosalia  shelf.  An  average  of  historic  hydrographic  data  is 
used  to  obtain  a  single  profile  of  stratification  below  200-m  depth  for  all  locations.  This  profile  is 
specified  as  a  decaying  exponential  with  a  length  scale  of  500  m.  Near-bottom  velocities  arc  not 
available  to  empirically  determine  the  bottom  friction  parameter  r  ( a is  proportional  to  r,  sec 
Clarice  and  Van  Gorder,  1986).  A  constant  value  of  0.05  cm/s,  used  in  previous  wave  studies 
(c.g.,  Chapman,  1987;  Church  et  al.,  1986),  is  chosen  for  all  sections. 

The  free  wave  parameters  computed  along  the  Mexican  coast  are  pictured  in  Fig.  4  for 
the  first  four  modes.  Mode  1  free  wave  phase  speeds  are  in  agreement  with  previous  values  for 
this  region  computed  by  EA,  except  near  Salina  Cruz.  The  close  agreement  between  our  mode  2 
and  EA’s  mode  1  suggests  that  they  may  have  missed  the  lowest  mode  in  their  analysis. 
Alongshore  variation  in  the  parameters  reflects  changes  in  shelf  width.  Wide  shelf  regions,  such 
as  near  Salina  Cruz  and  Mazatlan,  have  higher  free  wave  phase  speeds  while  narrow  shelf 
regions,  such  as  between  Puerto  Angel  and  Manzanillo,  have  weaker  wind  coupling  coefficients. 
Frictional  decay  length  scales,  measured  as  the  inverse  of  ajj ,  arc  shortest  in  wide  shelf  regions 
and  in  the  Gulf  of  California,  except  at  the  narrow  shelf  near  Santa  Rosalia.  The  free  wave 
parameters  are  reasonably  stable  to  changes  in  input  bathymetry  and  stratification  that  reflect  the 
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uncertainty  in  specifying  a  single  profile  for  a  given  section.  Mode  1  parameters  typically  do  not 
vary  by  more  than  10%  while  mode  4  parameters,  particularly  the  a ^ ,  can  change  by  50%  or 
more  in  some  cases.  Mode  4  is  the  highest  mode  considered  in  this  study.  As  noted  by  EA,  it  is 
important  to  specify  the  Baja  California  peninsula  for  the  Gulf  of  California  sections  rather  than 
assuming  a  deep  ocean  offshore.  The  presence  of  the  Baja  California  peninsula  increases  phase 
speeds  by  approximately  10-15%. 

The  free  wave  offshore  structures  for  pressure,  alongshelf  current,  and  density  at  section 
9  are  shown  in  Fig.  5.  The  profiles  reflect  the  hybrid  nature  of  free  waves  in  this  region  with  a 
more  barotropic  structure  over  the  shelf,  and  baroclinic  structure  over  the  slope.  The  offshore 
decay  scales,  as  noted  by  EA,  are  typically  less  than  1/2  the  width  of  the  gulf.  Using  a  flat  bottom 
offshore  condition  does  not  significantly  alter  the  wave  structures  over  the  mainland  shelf  and 
slope. 


4.  FORCED  WAVE  RESULTS 

Since  we  are  interested  in  storm  generated  events  rather  than  a  complete  hindcast  of  gulf 
shelf  observations,  it  is  desirable  to  statistically  isolate  a  propagating  signal  in  the  moored  data,  to 
compare  with  the  model  results.  CDG  and  EA  are  able  to  detect  poleward  propagation  in  coastal 
sea  level  records  from  Acapulco  to  the  Gulf  of  California  in  both  the  time  and  frequency 
domains.  We  prefer  to  isolate  the  propagating  signal  in  the  time  domain,  however,  since  wave 
events  are  non-dispersive,  short-term,  isolated  features  that  are  wide-banded  in  frequency  (0.02- 
0.37  cpd  from  EA).  One  seemingly  appropriate  time  domain  statistical  technique  is  complex 
empirical  orthogonal  function  (CEOF)  analysis  (Barnett,  1983).  CEOFs  have  been  used  to  iden¬ 
tify  propagating  variability  in  large  data  sets  (e.g.,  White  1987)  which  cannot  be  simply  accom¬ 
plished  with  a  standard  EOF  analysis  (Lorenz,  1956).  Tests  of  CEOF  analysis  on  synthetic  time 
scries  that  approximate  the  propagating  signals  of  this  study  indicate,  however,  that  the  method 
may  not  be  appropriate  for  wide-banded,  non-dispersive  processes  (see  chapter  3).  Instead,  we 
proceed  by  using  standard  EOFs  on  time  shifted  scries.  Coastal  adusted  sea  level  and  pressure 
observations  from  Acapulco  to  Isla  Tiburon  are  band-passed  to  include  only  frequencies  with 
poleward  phase  propagation  identified  by  EA,  and  shifted  in  time  to  obtain  maximum  correlation 
with  adjusted  sea  level  at  Guaymas  (Fig.  6).  An  EOF  of  these  shifted  time  series  has  a  lowest 
mode  which  describes  77%  of  the  variability  with  nearly  equal  spatial  amplitude  at  each  of  the 
four  locations.  When  mode  1  is  removed  from  each  time  scries,  residual  propagating  features  arc 
not  observed  (Fig  6).  Thus,  we  assume  that  the  propagating  signal  is  represented  by  this  single 
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mode.  The  propagation  speeds  derived  from  the  time  lags  shown  in  Fig.  6,  are  consistent  with  the 
previous  estimates  of  EA  obtained  from  sea  level  cross-spectra  (MW). 

Next,  a  multivariate  regression  analysis  is  performed  on  the  gulf  moored  observations  to 
estimate  propagating  and  locally  forced  variability.  Each  estimate  consists  of  two  inputs:  local 
forcing  represented  by  wind  stress  measured  at  W5  (equivalent  to  the  dominant  wind  EOF  mode 
described  in  MW),  and  the  remotely  forced  propagating  signal  described  above.  The  inputs  are 
lagged  in  time  to  coincide  with  maximum  correlation  with  the  estimand,  and  all  time  senes  are 
band-passed  filtered  (0.03-0.33  cpd).  Temperature  is  converted  to  density  using  the  relationship 
given  in  MW.  Density  at  Guaymas  has  slightly  higher  correlations  with  the  wave  signal  at  small, 
non-zero  time  lags.  In  this  study  the  zero-lagged  response  is  investigated. 

For  all  observations  in  the  gulf,  the  local  wind  mode  does  not  describe  a  significant  por¬ 
tion  of  the  variability.  This  is  true  for  lags  of  up  to  2  days  between  wind  stress  and  the  variable 
being  estimated.  MW  attribute  this  lack  of  local  wind  forcing  to  the  combination  of  weak  winds 
and  high  stratification  found  in  the  gulf  region.  In  contrast,  the  propagating  signal  accounts  for 
much  of  the  coastal  pressure  variability  along  the  coast  from  Salina  Cruz  to  Isla  Tiburon  (Table 
1).  The  amplitude  and  percentage  of  variance  explained  by  the  propagating  signal  is  weaker  at 
Salina  Cruz  than  at  mainland  coastal  stations  to  the  north,  consistent  with  wave  forcing  occurring 
north  of  Salina  Cruz  (EA).  The  offshore  pressure  structure  associated  with  the  propagating  signal 
decays  offshore  as  can  be  seen,  for  example,  at  Guaymas.  The  propagating  amplitude  at  Topolo- 
bampo  is  lower  than  that  at  neighboring  stations  due  to  the  offshore  location  (100-m  isobath)  of 
the  sensor.  The  propagating  signal  is  also  detected  in  pressure  observations  on  the  Baja  Califor¬ 
nia  shelf,  but  with  low  amplitude.  The  signal  on  this  side  of  the  gulf  lags  the  mainland  coast  sig¬ 
nal  by  1.3-1. 8  days  at  Santa  Rosalia-Guaymas,  and  2.4  days  at  San  Francisquito-Isla  Tiburon. 

In  the  Gulf  of  California,  density  and  alongshelf  current  measurements  also  have  a 
detectable  propagating  component.  Density  fluctuations  associated  with  the  propagating  signal 
are  significant  below  the  upper  mixed  layer  at  Guaymas  and  Topolobampo  and  explain  approxi¬ 
mately  40  to  60%  of  the  variability  (Table  1).  Density  fluctuations  at  and  below  70-m  depth  at 
Santa  Rosalia  are  also  correlated  with  the  propagating  signal  at  an  approximate  9  day  lag  relative 
to  Guaymas,  as  noted  in  MW.  The  propagating  signal  explains  a  significant,  although  low  per¬ 
centage,  of  the  alongshelf  current  variability  on  the  mainland  coast.  MW  report  that  the  flow 
field  is  more  complex  than  the  pressure  and  density  fields  with  considerable  variability  not 
directly  accounted  for  by  either  local  winds  or  propagating  events.  Thus  individual  current  time 
scries  will  be  directly  compared  with  model  hindcasts. 
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Since  observed  propagation  speeds  are  most  similar  to  mode  1  free  wave  phase  speeds 
(EA,  CDG),  the  empirically  derived  propagating  signal  is  first  compared  with  mode  1  model 
results,  followed  by  an  examination  of  the  effects  of  including  higher  modes.  The  comparison 
between  observed  and  mode  1  coastal  pressure  is  shown  in  Fig.  7.  In  this  initial  run,  is 

set  to  zero  under  the  hypothesis  that  most  of  the  wave  variability  is  generated  north  of  Salina 
Cruz.  The  model  does  reasonably  well  in  reproducing  the  largest  1984  events  from  Topolo- 
bampo  to  Isla  Tiburon  which  are  associated  with  hurricanes  Fausto  and  Genevieve  in  July  (Fig. 
8).  Two  of  the  three  smaller  events  in  late  July  and  August  associated  with  weaker  tropical 
storms  are  also  reproduced  by  the  model,  although  the  amplitudes  are  underestimated  by  approxi¬ 
mately  30-50%.  The  model  does  predict  an  event  in  early  June  that  is  not  apparent  in  the  obser¬ 
vations.  Maximum  correlations  between  the  EOF  propagating  signal  and  the  mode  1  coastal 
pressure  hindcasts  indicate  that  observed  arrival  time  on  the  mainland  shelf  is  36  hours  earlier 
than  the  model  estimate  at  Toplobampo,  24  hours  at  Guaymas,  and  18  hours  at  Isla  Tiburon. 
Correlations  between  the  EOF  signal  and  the  mode  1  hindcast  are  approximately  0.6  in  the  Gulf 
of  California. 

The  model  indicates  that  wave  generation  should  occur  primarily  between  Manzanillo 
and  Mazatlan  where  storms  are  strongest  and  closest  to  shore  (Fig.  8).  The  corresponding  events 
in  observed  sea  level,  however,  appear  much  further  south.  For  example,  the  strong  July  events 
are  present  only  in  the  model’s  mode  1  coastal  pressure  north  of  Manzanillo  in  contrast  to  the 
observations  which  seem  to  show  the  events  as  far  south  as  Salina  Cruz  (see  also  EA  for  similar 
events  in  1971,  1973-1974).  Storm  tracks  do  not  show  either  Fausto  or  Genevieve  far  enough 
south  to  generate  the  event  at  Salina  Cruz  using  our  extrapolation  technique.  The  sea  level  events 
associated  with  weaker  tropical  storms  in  late  July  and  August  also  appear  further  south  than  the 
model  indicates.  The  model  fit  is  not  improved  by  using  Salina  Cruz  adjusted  sea  level  to  specify 
<t>i(0,  r).  For  the  value  of  bottom  friction  used  in  this  study,  the  signal  input  at  Salina  Cruz 
decays  by  90%  before  reaching  Acapulco,  while  the  observations  show  little  decay  over  this  dis¬ 
tance.  It  is  not,  however,  simply  a  matter  of  adjusting  the  friction  parameter  as  the  free  wave 
speeds  are  too  slow  for  the  event  seen  at  Salina  Cruz  to  reach  Acapulco  in  the  observed  time 
span.  The  propagation  speed  between  Salina  Cruz  and  Acapulco  during  the  large  July  events  is 
almost  8  m/s,  nearly  three  times  as  large  as  the  mode  1  phase  speed.  This  observed  speed,  how¬ 
ever,  is  still  far  too  slow  for  a  barotropic  Kelvin  wave  which  has  a  speed  on  the  order  of  200  m/s. 
This  discrepancy  in  the  location  of  wave  generation  is  most  likely  due  to  the  inherent  limitations 
of  estimating  the  coastal  wind  field.  The  radial  profiles  of  tangential  winds  arc  for  wcll- 
dcvclopcd  storms.  During  the  early  formation  stages  of  these  storms,  which  often  occur  in  the 
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vicinity  or  south  of  Salina  Cruz,  the  wind  field  is  certainly  far  less  structured  and  may  not  be  ade¬ 
quately  represented  by  our  simple  coastal  wind  extrapolation.  We  could  not  find,  however,  a  sim¬ 
ple  variation  of  the  radial  storm  profile  that  is  more  consistent  with  all  the  observed  events. 
Nonetheless,  the  temporal  changes  of  propagating  events  in  observed  sea  level  between  Acapulco 
and  the  Gulf  of  California  (Fig.  7)  do  seem  to  support  the  model  result  that  some  forcing  must 
occur  north  of  Acapulco.  In  addition,  the  model  comparisons  north  of  Acapulco  verify  that  the 
typical  winds  around  storms  in  this  region  are  sufficient  to  generate  the  observed  sea  level  ampli¬ 
tudes.  The  event  of  roughly  15  cm  in  July  is  reproduced  by  the  model  when  forced  by  coastal 
winds  of  20-30  m/s  near  Manzanillo  and  Mazatlan.  When  PAMII  winds  with  rms  speeds  of  5  m/s 
are  used  to  force  the  model  in  the  gulf,  the  associated  response  is  on  the  order  of  1  cm  for  sea 
level  and  1  cm/s  in  current  speed  again  confirming  the  weak  local  wind  forcing  in  the  gulf. 

The  model’s  mode  1  pressure  and  density  cross-shelf  structures  at  Guaymas  are  qualita¬ 
tively  similar  with  profiles  obtained  from  the  multivariate  analysis  (Fig.  9).  In  addition  to  the 
measured  100-m  and  200-m  bottom  pressure  and  adjusted  coastal  sea  level,  pressure  is  computed 
at  various  depths  above  the  200-m  pressure  sensor  using  density,  obtained  from  the  thermistor 
chain,  and  the  hydrostatic  equation  (described  in  MW).  The  mode  1  and  observed  pressure  and 
density  fields  are  normalized  to  correspond  to  a  10  mbar  pressure  rise,  or  10  cm  sea  level  rise,  at 
the  coast.  The  two  pressure  fields  have  largest  amplitudes  at  the  coast,  decaying  offshore  and 
with  depth.  The  observed  pressure  field,  however,  as  a  stronger  baroclinic  component  and  decays 
approximately  twice  as  fast  offshore.  The  mode  1  and  observed  cross-shelf  density  fields  are  also 
qualitatively  similar  in  structure  with  largest  fluctuations  in  the  thermocline  region  between  40-m 
and  70-m  depth  with  amplitudes  stronger  near  the  shelf  break  than  at  the  coast  (Fig.  10).  The 
model’s  mode  1  amplitudes,  however,  are  almost  3  times  weaker  than  those  observed.  As  model 
pressure  and  density  are  in  hydrostatic  balance,  the  barotropic  structure  of  the  mode  1  pressure 
field  is  a  direct  result  of  the  model’s  weak  density  response  over  the  shelf. 

On  the  Santa  Rosalia  shelf,  the  density,  alongshelf  current,  and  pressure  fluctuations  are 
found  to  be  related  to  the  propagating  signal  (Table  2),  but  at  different  lags.  The  propagating  sig¬ 
nal  in  pressure  is  lagged  roughly  1.5  days  relative  to  Guaymas  pressure,  while  density,  at  70-m 
and  100-m  depth  at  the  100-m  isobath,  is  lagged  by  8.5  days.  The  mode  1  travel  time  around  the 
head  of  the  gulf  is  roughly  8.2  days  with  the  next  highest  mode  taking  15.5  days,  indicating  that 
wave  energy  exhibits  mode  1  propagation  even  over  the  very  complicated  bathymetry  of  the 
northern  gulf.  The  model’s  mode  1,  however,  has  the  strongest  density  fluctuations  in  the  ther¬ 
mocline  region  above  40-m  depth  in  contrast  to  the  observations  which  only  show  a  detectable 
propagating  response  below  70  m.  The  model’s  mode  1  also  underestimates  the  amplitude  of  the 


69 


deep  density  fluctuation  by  two  orders  of  magnitude.  For  mode  1  to  reproduce  the  observed  den¬ 
sity  amplitude  would  require  an  unrealistic  50-cm  rise  in  coastal  sea  level.  The  time  lag  of  the 
Santa  Rosalia  pressure/current  fluctuations  relative  to  Guaymas  is  much  too  short  for  this 
response  to  be  associated  with  the  equatorward  signal  observed  in  density.  The  pressure/currcnt 
signal  at  Santa  Rosalia  would  appear  to  be  the  offshore  edge  of  the  wave  as  it  travels  poleward 
along  the  mainland  shelf.  The  decay  scale  across  the  gulf  (150  km),  however,  does  not  agree 
with  the  offshore  decay  seen  on  the  Guaymas  shelf  (30  km)  or  with  the  model’s  mode  1  pressure 
(50  km).  In  addition,  the  time  lag  between  Guaymas  and  Santa  Rosalia  pressure  cannot  be 
accounted  for  by  the  first  wave  mode.  That  the  Santa  Rosalia  pressure/current  signal  is  not  sim¬ 
ply  an  extention  of  the  wave  across  the  gulf  is  also  supported  by  the  insignificant  correlation 
between  Santa  Rosalia  alongshelf  currents  and  Guaymas  sea  level  during  the  1983  summer 
deployment.  In  constrast,  Santa  Rosalia  density,  which  we  do  directly  associate  with  the  wave 
events,  is  correlated  in  both  summers.  MW  suggest  that  the  current  signal  at  Santa  Rosalia  may 
be  associated  with  eddy-like  features  setup  by  the  larger  wave  events  rather  than  a  component  of 
the  wave  structure. 

Alongshore  current  observations  on  the  mainland  shelf  agree  with  the  model’s  mode  1 
vertical  structure,  but  not  in  cross-shelf  structure.  During  the  energetic  July  wave  events, 
alongshelf  currents  at  Topolobampo  are  in  good  agreement  with  the  magnitude  and  direction  of 
the  model’s  mode  1  response  and  confirm  that  these  currents  are  essentially  barotropic  over  the 
shelf  (Fig.  1 1).  The  comparison  of  the  weaker  currents  in  August,  however,  is  poorer  but  may  be 
due  to  the  masking  of  wave  variability  by  other  processes.  This  appears  to  be  the  case  at  Guay¬ 
mas  (Fig.  12)  where  considerably  more  variability  occurs  which  is  not  associated  with  storm 
events.  The  July  wave  event  is  assumed  strong  enough,  however,  to  allow  a  direct  examination 
of  currents  during  this  time  period.  Currents  measured  in  total  water  depths  of  200  m,  50  m,  and 
15  m  show  a  clear  diminishing  of  current  speed  towards  the  coast,  in  contrast  to  mode  1  currents 
which  increase  towards  the  coast.  Observed  currents  at  M4  (15-m  depth)  have  peak  speeds  of 
10  m/s  compared  to  model  predictions  of  over  50  m/s.  In  addition  the  direction  of  15-m  currents 
appears  to  be  counter  to  observed  offshore  currents  and  model  currents.  There  is  also  a  slight 
phase  lag  between  peak  current  occurrences  across  the  shelf  during  the  July  propagating  events, 
M6  leads  M8  by  18  hours,  that  is  not  explained  by  mode  1.  Although  the  exact  details  of  the 
complicated  Guaymas  flow  field  may  be  beyond  the  scope  of  this  wave  model,  the  model’s 
increase  in  current  amplitude  towards  the  coast  is  certainly  unrealistic. 

While  mode  1  is  successful  in  accounting  for  most  of  the  propagating  events  seen  in 
coastal  sea  level  in  the  gulf,  the  cross-shelf  structures  of  alongshelf  current,  density,  and  pressure 
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suggest  the  presence  of  higher  modes.  Specifically,  the  model’s  alongshelf  current  field  should 
decrease  rather  than  increase  towards  the  coast,  and  the  amplitude  of  density  fluctuations  are 
underpredicted  by  a  factor  of  three,  thus  underpredicting  the  baroclinic  component  to  the  pressure 
field  over  the  shelf.  Several  model  runs  are  made  to  investigate  the  effect  of  including  higher 
modes  in  the  model  hindcast. 

The  first  set  of  model  runs  use  all  four  modes  shown  in  Fig.  5,  without  frictional  cou¬ 
pling  between  modes  ( a ij  =  0  for  i  *  j ).  The  energy  in  modes  3  and  4  is  essentially  dissipated 
before  reaching  the  Gulf  of  California  for  the  chosen  magnitude  of  the  frictional  parameter. 
Mode  2,  with  a  slower  propagation  speed,  arrives  in  the  gulf  approximately  5  days  after  mode  1. 
This  result  is  inconsistent,  however,  with  the  inability  to  detect  higher  propagating  signals  in  coa¬ 
stal  sea  level  observations  (Fig.  6).  A  similar  result  is  found  in  the  fully  coupled  case.  Thus,  we 
assume  that  higher  modes  (j  >  1)  do  not  freely  propagate  into  the  gulf,  but  are  excited  by  the 
dominant  mode  1  through  the  frictional  coupling  as  it  travels  along  the  coast.  In  the  next  set  of 
model  runs,  higher  modes  are  excited  by  mode  1  only  at  section  9  where  the  observations  arc 
made.  For  2  coupled  modes  at  Guaymas,  the  velocity  field  is  improved  somewhat  with  a  less 
rapid  increase  in  current  amplitude  towards  the  coast.  Density  fluctuations,  however,  are  now 
even  weaker  over  the  shelf.  Modes  3  and  4  are  only  weakly  coupled  with  mode  1  at  Guaymas 
and  do  not  alter  the  cross-shelf  structures  obtained  from  the  first  two  modes. 

These  problems  can  be  addressed  by  increasing  the  frictional  parameter  which  would 
damp  out  modes  higher  than  mode  1  from  freely  propagating  into  the  gulf,  as  indicated  by  the 
observations,  and  increase  mode  coupling  which  might  explain  the  discrepancies  between  mode  1 
and  the  observations.  This  increase  in  friction,  however,  would  also  dissipate  energy  in  mode  1, 
which  is  contrary  to  the  minimal  decay  in  wave  amplitude  between  Guaymas  and  Isla  Tiburon 
(Table  1),  and  would  result  in  even  less  available  energy  to  increase  the  amplitude  of  density 
fluctuations. 


5.  DISCUSSION  AND  SUMMARY 

Propagating  sea  level  events  in  the  Gulf  of  California  are  similar  to  mode  1  coastal- 
trapped  waves.  The  hindcast  of  mode  1  coastal  pressure  using  a  simple  estimate  of  storm  winds 
at  the  coast,  is  successful  in  reproducing  the  amplitude  and  approximate  propagation  speed  of 
energetic  wave  events  north  of  Mazatlan.  Wave  generation,  however,  is  predicted  to  occur  pri¬ 
marily  between  Manzanillo  and  Mazatlan  during  1984,  while  observations  indicate  generation 
occurring  at  least  as  far  south  as  Salina  Cruz  for  some  events.  This  discrepancy  is  most  likely 
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due  to  our  limited  knowledge  of  the  coastal  wind  field  rather  than  an  inherent  problem  with  the 
model.  At  the  Guaymas  shelf,  the  observed  pressure  field  is  trapped  to  the  coast  in  the  same 
manner  as  mode  1,  although  the  observed  field  is  more  baroclinic  over  the  shelf  and  has  a  shorter 
offshore  decay  scale.  The  observed  and  mode  1  density  fields  at  Guaymas  are  also  similar  in 
structure,  although  observed  amplitudes  are  nearly  three  times  larger  than  predicted  by  mode  1. 
Alongshelf  currents  are  barotropic  during  energetic  wave  events  at  the  midshelf,  in  accordance 
with  mode  1.  Currents  are  observed  to  decay  in  amplitude  towards  the  coast,  however,  in  contrast 
to  mode  1  current  profiles  which  have  highest  amplitudes  at  the  coast.  Higher  modes  reduce  the 
increase  in  current  speed  towards  the  coast,  but  at  the  expense  of  less  realistic  pressure  and  den¬ 
sity  profiles  over  the  shelf.  Moreover,  alongshelf  current  observations,  made  primarily  in  depths 
shallower  than  20  m,  are  probably  in  frictional  regions  unsuitable  for  the  coastal-trapped  wave 
model  dynamics.  Thus  the  use  of  higher  modes  to  hindcast  currents  outside  the  model  domain  is 
not  justifiable.  The  inclusion  of  modes  higher  than  mode  4  are  unlikely  to  significantly  improve 
the  hindcast  results,  while  increasing  the  frictional  parameter  to  increase  mode  coupling,  would 
also  increase  the  dissipation  of  the  wave  events  in  contrast  to  the  observations. 

We  have  demonstrated  that  coastal  winds  during  a  hurricane  or  tropical  stom  in  the 
eastern  north  Pacific  are  sufficient  to  excite  observed  wave  amplitudes.  How  far  south  the  gen¬ 
eration  occurs  is,  however,  unresolved  as  it  is  uncertain  what  the  nature  of  the  signal  is  at  Salina 
Cruz  where  storms  are  generally  weak.  CDG  and  EA  report  low  correlations  or  coherences 
between  sea  level  at  Salina  Cruz  and  more  northern  stations  for  the  summers  of  1971  and  1973- 
1974,  which  supports  their  argument  that  forcing  occurs  north  of  Salina  Cruz.  While  propagating 
events  are  weak  or  absent  at  Salina  Cruz  in  the  1971  summer,  several  energetic  events  do  seem  to 
appear  as  far  south  as  San  Jose  (shown  in  Fig.  3)  in  1973-1974  (see  EA  Figs.  2  and  4).  A  better 
description  of  the  wind  field  along  the  Mexican  coast  is  needed,  but  particularly  south  of  Salina 
Cruz,  in  order  to  uncquivocably  link  all  propagating  events  to  tropical  storms  and  hurricanes. 

In  addition,  the  results  from  this  study  suggest  that  the  coastal-trapped  wave  model  docs 
not  accurately  estimate  the  amplitude  of  the  density  response.  Chapman  (1987)  and  Brink  (1982) 
also  report  poor  hindcasts  of  density,  but  contend  that  mixed  layer  and  other  non-wave  ;  'latcd 
phenomena  contaminate  the  wave  signal  in  density.  In  this  study,  we  are  able  to  cmpincally 
identify  a  propagating  signal  in  density  independently  of  the  model  and  find  that  lower,  more 
dominant  modes  are  not  able  to  explain  the  observed  response.  The  model  dynamics  assume  that 
density  fluctuations  are  due  solely  to  the  vertical  advection  of  the  mean  stratification.  This 
assumption  is  consistent  with  the  description  of  the  temperature  field  in  the  gulf  presented  by 
MW.  The  mean  stratification  does  vary  over  the  1984  instrument  deployment  in  contradiction  to 
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a  model  assumption,  but  the  magnitude  is  not  significant  enough  to  explain  the  model’s  underes¬ 
timate  of  density.  The  problem  seems  to  involve  the  model’s  underestimate  of  the  vertical  velo¬ 
city.  If  so,  the  addition  of  more  realistic  surface  and  bottom  boundary  layers  may  improve  the 
performance  of  the  model. 
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TABLE  1.  Results  from  the  multiveriate  regression  of  the  empirically  determined  propagating  signal  and  alongshelf  wind  stress 
on  moored  observations  in  the  Gulf.  The  alongshelf  wind  stress  does  not  have  a  significant  component  in  the  observations  and  is 
not  listed. 


Instrument  Location 

a 

Lag 

(days) 

% 

Variance 

PRESSURE 

Salina  Cruz  (SC) 

0.70  ±0.11 

-8.5 

31.7* 

Acapulco  (AC) 

0.97±0.25 

-6.4 

61.9* 

Topolobampo  100  m  (TO) 

0.63  ±  0.04 

-1.2 

78.9* 

Guaymas  Harbor  (GY) 

1.20  ±0.04 

0.0 

85.4* 

Guaymas  100  m  (P4) 

0.5910.02 

0.0 

67.7* 

Guaymas  200  m  (P5) 

0.25  ±  0.06 

0.0 

25.9* 

Isla  Tiburon  7  m  (IT) 

1.0310.02 

0.9 

83.3* 

S.  Francisquito  7  m  (SF) 

0.38 1 0.04 

3.3 

19.8* 

S.  Rosalia  Harbor  (SR) 

0.29 1 0.02 

1.3 

36.6* 

S.  Rosalia  100  m  (P6) 

0.27 1 0.08 

1.8 

20.5* 

DENSITY 

Ml  70/100  m 

-0.37 1 0.06 

-1.1 

24.4* 

M4  2/15  m 

0.05 1 0.01 

0.0 

3.6 

14/15  m 

-0.06 1 0.04 

0.0 

1.7 

M6  10/50  m 

0.02 1 0.02 

0.0 

0.4 

20/50  m 

-0.1110.04 

0.0 

4.6 

30/50  m 

-0.38 1 0.08 

0.0 

18.7* 

M7  100/100  m 

-0.54 1 0.02 

0.0 

53.7* 

M8  10/200  m 

-0.0310.06 

0.0 

0.7 

20/200  m 

-0.2010.04 

0.0 

8.8* 

30/200  m 

-0.4810.10 

0.0 

19.2* 

40/200  m 

-0.7410.13 

0.0 

30.5* 

50/200  m 

-0.9110.11 

0.0 

45.0* 

60/200  m 

-0.9210.06 

0.0 

58.7* 

70/200  m 

-0.8410.04 

0.0 

62.5* 

85/200  m 

-0.6510.08 

0.0 

58.2* 

1 15/200  m 

-0.4010.02 

0.0 

56.9* 

1 30/200  m 

-0.3010.02 

0.0 

53.8* 

145/200  m 

-0.2410.02 

0.0 

51.3* 

1 60/200  m 

-0.2110.01 

0.0 

49.0* 

175/200  m 

-0.1810.01 

0.0 

49.8* 

200/200  m 

-0.11  ±0.02 

0.0 

47.4* 

M9  50/1000  m 

0.1010.02 

0.0 

2.2 

Mil  10/100  m 

0.05 1 0.08 

0.0 

0.8 

30/100  m 

0.0710.15 

0.0 

0.9 

70/100  m 

-0.2410.04 

9.6 

31.3* 

100/100  m 

-0.25  1 0.02 

8.5 

38.3* 

ALONGSHELF  CURRENT 

Ml  70/100  m 

0.90 1 0.23 

-1.6 

8.3* 

M4  5/15  m 

0.1510.02 

-1.8 

8.4* 

M6  20/50  m 

-0.9510.12 

2.1 

12.0* 

M8  10/200  m 

2.7110.54 

-0.25 

16.2* 

M9  50/1000  m 

1.0610.42 

2.1 

7.9 

Mil  10/100  m 

-1.1610.25 

1.7 

16.6* 

30/100  m 

-1.2410.17 

1.5 

25.0* 

70/100  m 

-I.l510.12 

1.5 

21.6* 

*  above  95%  significance  level 

The  regression  coefficient,  a  (g/cm3/mh  for  density,  cm /slmb  for  current),  is  relative  to  the  magnitude  of  the  propagating  signal  at 
Guaymas.  Time  series  are  shifted  in  time  (lag)  to  obtain  maximum  correlation  with  the  propagating  signal.  Pressure  observations 
which  do  not  include  a  depth  refer  to  adjusted  coastal  sea  level. 


Figure  I. 


The  Gulf  of  California  showing  mooring  (M),  bottom  pressure  (P),  and  PAMII  meteoro 
logical  station  (W)  locations. 
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w(0  "  wmoxexP[-(r-rmax)/A] 
rmox  "  30  n-m- 

A  =  630exp[-^mox  -  80)/50] 

Figure  2. 

The  extrapolation  of  tropical  storm  and  hurricane  winds  to  the  coasts,  based  on  the  radial 
wind  speed  profiles  of  Riehl,  1954.  The  exponential  decay  scale  for  hurricane  winds  is  a 
function  of  wind  speed.  The  wind  field  is  taken  to  decay  linearly  for  tropical  storms. 
w=wind  speed,  w^^maximum  wind  speed,  i^radial  distance  from  storm  center, 
r  =radius  of  maximum  wind  speed. 


78 


Figure  3. 

The  bathymetry  of  the  Pacific  coast  of  Mexico  showing  the  200-m  (dashed)  and  1000-, 
2000- ,  and  3000-m  isobaths  (solid).  Inset  is  a  map  showing  the  12  sections  for  the 
coastal- trapped  wave  model.  Original  figure  from  Enfield  and  Allen  (1983). 
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Figure  6. 


(a)  Coastal  bottom  pressure  at  Isla  Tiburon  (7-m  depth),  and  Topolobampo  (100-m),  and 
adjusted  sea  level  at  Guaymas  and  Acapulco  with  each  series  time  shifted  to  obtain  max¬ 
imum  correlation  with  the  Guaymas  time  series.  The  dominant  EOF  mode  of  the  four  time 
series  is  plotted  last,  (b)  The  same  four  time  series  minus  mode  1.  ine  residual  time 
series  are  not  significantly  correlated  and  do  not  exhibit  propagating  variability. 


1984  OBSERVED  COASTAL  PRESSURE 
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Tracks  of  Eastern  North  Pacific  tropical  storms  for  1984  from  Gunther  and  Cross  (1985). 
All  storms  reached  hurricane  intensity  along  portions  of  the  coast  except  Julio. 


OFFSHORE  (KM) 


OFFSHORE  (KM) 


Figure  9. 

The  Guaymas  cross-shelf  structure  of  pressure,  in  millibars,  from  the  multi-variate  regres¬ 
sion  analysis  (observed)  and  the  coastal-trapped  wave  model  mode  1.  The  pressure  fields 
are  normalized  to  10  millibars  at  the  coast. 
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Figure  10. 


The  Guaymas  cross-shelf  structure  of  density,  in  g/cm3,  from  the  multi-variate  regression 
analysis  (observed)  and  the  coastal-trapped  wave  model  mode  1.  The  density  fields  are 
normalized  to  correspond  to  a  10  millibar  pressure  rise  at  the  coast. 


CHAPTER  3 


COMPLEX  EMPIRICAL  ORTHOGONAL  FUNCTION  ANALYSIS: 
WIDE  VERSUS  NARROW  BANDED  PROCESSES 


1.  INTRODUCTION 

Empirical  Orthogonal  Functions  (EOFs,  see  Lorenz,  1956  for  a  description)  arc  now  a 
standard  technique  for  describing  coherent  variability  in  geophysical  data  sets.  EOFs  arc 
designed  to  give  the  most  efficient  description  of  the  observed  variability,  often  reducing  large 
data  sets  to  a  few  dominant  modes.  These  dominant  modes  are  sometimes  useful  in  representing 
physical  process  in  the  data  (e.g.,  Kundu  et  al.,  1975),  although  the  requirement  of  orthogonality 
in  both  time  and  space,  the  use  of  the  covariance  (or  correlation)  matrix  at  zero  lag  to  describe  the 
variability,  and  the  particular  sampling  of  the  observed  field  require  careful  consideration  before 
assigning  any  physical  significance  to  a  given  mode.  In  particular,  the  use  of  the  zero-lagged 
covariance  matrix  limits  the  analysis  to  the  identification  of  standing  wave  patterns. 

An  alternative  method  for  detecting  propagating  patterns  is  complex  empirical  orthogo¬ 
nal  function  analysis  (CECFs,  Rasmusson  et  al.,  1983).  This  method,  carried  out  in  the  time 
domain,  has  been  used  to  detect  large-scale  travelling  patterns  in  meteorological  (Barnett,  1983) 
and  oceanographic  observations  (White  and  Tabata,  1987).  Horcl  (1984)  discusses  several  tests 
of  the  method  on  simulated  data.  In  this  chapter,  we  investigate  the  suitability  of  CEOFs  for  iso¬ 
lating  propagating  signals  associated  with  coastal-trapped  waves  along  the  Mexican  Pacific  coast. 
Specifically,  we  are  interested  in  isolating  a  single  CEOF  mode  which  characterizes  the  propaga¬ 
tion  observed  in  coastal  sea  level  and  moored  observations  in  the  Gulf  of  California.  After 
several  test  cases  using  simulated  data,  we  conclude  that  the  method  is  not  appropriate  for  this 
particular  application,  and  may  not  be  suitable  in  general  for  broad-banded  processes.  A  brief 
description  of  the  method  is  given  in  section  2,  the  model  test  is  described  in  section  3,  and  a 
summary  follows  in  section  4. 
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2.  METHOD 

CEOFs  involve  obtaining  quadrature  information  about  a  scalar  field  without  resorting 
to  a  frequency  domain  analysis.  Following  Barnett  (1983),  we  represent  a  scalar  field 
where  j  signifies  spatial  position  (j  =  1,2,...,  N)  and  t  is  time,  in  the  usual  Fourier  sense  as 

Ujit)  =  (oo)cos  (03  r)  +  (co)sin(co  r )] .  (1) 

CO 

To  obtain  quadrature  information  between  stations,  a  complex  representation  of  Uj  (t )  is  invoked, 

Uj(t)  =  ZCj(C£>)eia>‘  ,  (2) 

CO 

where  Cj  (co)  =  (gj)  +  ibj  (co).  Then  Uj  it)  can  be  written  in  the  form, 

Ujit)  =  X  {[ttj  (o3)cos(cor )  +  &y(co)sin(cor)] 

CO 

+  i  [bj  (co)cos(cof )  -  cij  (co)sin(cor )]} 

=  Uj(t)  +  iUj(t)  (3) 

The  real  part  of  (3)  is  the  original  scalar  field;  the  imaginary  part  is  the  quadrature  function  of 
Ujit),  or  the  Hilbert  transform  (Thomas,  1969).  As  is  evident  from  the  Fourier  representation, 
the  Hilbert  transform  phase  shifts  each  spectral  component  of  Uj  it)  by  90°  in  time.  In  this  study, 
u  is  obtained  from  the  Fourier  coefficients. 

The  covariance  matrix  of  Uj  it)  is 

C jk  =  <U j  *  it  )Ukit  )>t  (4) 

where  <...>,  denotes  time  averaging.  It  follows  from  standard  spectral  analysis  theory  that 

co„ 

Cjk  =  j  <bjkia)du 


(5) 
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where  §jk  is  the  cross  spectrum  and  C0„  is  the  Nyquist  frequency.  Thus,  C jk  is  equivalent  to  the 
cross-spectral  matrix  integrated  over  all  frequencies. 

As  in  a  standard  EOF  analysis,  complex  eigenvectors,  Bn  ( x ),  are  obtained  from  C .  The 
complex  temporal  expansions  are  given  by, 

AnU)=XUj(t)Bn(x).  (6) 

7=1 

U  can  then  be  represented  as, 

Uj(t)=  ZAn(t)Bn*(x).  (7) 

n= 1 


As  an  example,  CEOFs  are  obtained  for  a  traveling  sinusoid  of  the  form 
u(x,t)  =  a  sin(£0;t -co0 1 ).  The  spatial  phase  is  defined  by 

0„  (x )  =  arctan  [ImBn  (x )  /  ReBn  ( x )]  .  (8) 

In  this  case,  there  is  only  one  eigenfunction  with  d01  /  dx  ~k0.  The  spatial  amplitude  is  given 
by 

5rtU)  =  [5„(x)fin*(x)]as.  (9) 

In  the  case  of  a  constant  amplitude  sinusoid,  S-l(x)=a.  The  temporal  phase  is  defined  by 

$n(t)  =  arctan[ImAn(t) / ReAn(t)]  .  (10) 

Again,  all  energy  is  in  the  first  mode  with  dQ^t)/  dt  =(£>0.  The  temporal  amplitude  is  defined 
by 

i?n(r)  =  [A„(OA„*(0]a5  (ID 

For  a  sinosoid  whose  amplitude  docs  not  vary  with  time,  R\(t)  is  a  constant.  It  follows  that  the 
phase  speed  is  obtained  by 


cn(x,t)  =  a/k  =(d$n(t)/dt)  /  (dQn(x)/dx) , 


(12) 
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which  in  the  case  of  the  sinusoid  is  c  j  (x ,  t )  =  C0o  I  ka. 


3.  TEST  CASE 

The  coastal-trapped  wave  test  case  (Fig.  1)  consists  of  a  series  of  irregularly  spaced 
Gaussian  bumps  that  propagate  non-dispersively  past  a  one-dimensional,  ten-station  array.  The 
time  series  schematically  represent  coastal  sea  level  observations  during  one  tropical  storm  sea¬ 
son  (6  months),  at  equally  spaced  stations  over  a  distance  equivalent  to  the  separation  between 
Salina  Cruz  and  Guaymas  (2500  km,  see  chapter  2,  Figure  3).  Wave  events  have  10  day  period, 
constant  amplitude,  and  2.5  m/s  phase  speed  (215  km/day).  As  in  the  case  of  the  travelling 
sinusoid  above  and  in  the  absence  of  any  background  noise,  we  seek  a  single  mode  which  charac¬ 
terizes  this  simple  travelling  pattern. 

CEOFs,  however,  are  unable  to  isolate  all  the  variability  in  a  single  mode.  The  first 
mode,  describing  72%  of  the  total  variance,  is  used  to  reconstruct  the  time  scries  at  each  station 
(Fig.  2).  The  wave  field  is  only  partially  accounted  for  by  mode  1,  with  the  percentage  of  vari¬ 
ance  described  at  each  station  decreasing  towards  the  ends  of  the  array  (Fig.  3).  The  average  tem¬ 
poral  phase  is  used  to  obtain  a  representative  to,  and  k  is  obtained  from  the  slope  of  the  spatial 
phase  function.  The  phase  speed  of  mode  1,  computed  from  equation  12,  underestimates  the  true 
phase  speed  by  33%. 

The  limitations  of  CEOFs  found  in  this  particular  test  case  arc  due  in  part  to  the  wide- 
banded  nature  of  the  wave  signal.  This  is  illustrated  in  the  following  example.  The  complex 
eigenvector  B  ,(.x)  obtained  from  the  covariance  matrix  C  has  an  arbitrary  phase.  The  absolute 
phase  is  determined  when  B  jOt )  is  projected  on  the  data  to  obtain  A  i(/)  from  equation  6.  Con¬ 
sider  a  spatial  phase  shift  of  B  j( x )  such  that  at  station  x$(j=5),  ©1(^5)  =  0.  Then  from  equa¬ 
tions  3  and  7, 

U  (x5,t )  =  u  (x 5,t )  +  iu  (. x5,t ) 


=  a  1  (t)b  1  (jx5)  +  i  [a\  ( t  )b \  (.x5)]  (13) 

where  a',  aln,  b„,  bln,  denote  the  real  and  imaginary  components  of  An  and  Bn,  and  b\  (x5)  =  0 
for  this  particular  phase  orientation  of  B  {(x ).  Hence,  the  real  part  of  A  j(r)  is  proportional  to  the 
observed  wave  signal  at  station  5,  and  the  imaginary  part  of  A  j(r)  is  proportional  to  the  Hilbert 
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transform  of  the  wave  signal.  The  real  and  imaginary  parts  of  A  ]  arc  pictured  in  Fig.  4.  At  the 
remaining  stations,  where  b\  (x,)  does  not  necessarily  equal  zero,  a  simple  linear  superposition 
of  a  j  and  a\  cannot  reproduce  liL  at  each  station.  In  the  case  of  the  sinusoid,  the  imaginary  pan 
of  A  i  is  simply  a  negative  cosine  function.  Thus,  the  real  and  imaginary  components  of  A  j,  with 
the  appropriate  spatial  phase  weighting,  can  recover  the  signal  entirely  at  any  station.  It  becomes 
less  likely  that  a  superposition  of  two  temporal  functions  can  reproduce  a  propagating  signal  as 
the  energy  in  that  signal  is  spread  out  more  in  frequency.  This  result  for  wide-banded  functions 
was  implied  by  Horcl  (1984)  for  a  similar  test  case.  He  notes  that  wide-banded  signals  arc  better 
represented  by  the  CEOF  as  the  phase  speed  increases.  In  this  limit,  the  temporal  lag  between 
time  scries  is  reduced  and  the  simple  superposition  discussed  above  is  more  tenable. 

An  additional  consideration  is  the  ability  of  CEOFs  to  reproduce  the  phase  speed,  c . 
For  a  non-dispersive  process,  c  is  a  constant  and  k  varies  linearly  with  frequency.  The  CEOF, 
however,  determines  a  single  (CO ,k)  pair.  The  spectrum  for  the  Gaussian  bump  wave  train  (Fig. 
5)  is  wide-banded  at  low  frequencies  (<  0.2  cpd),  making  a  determination  of  a  single  (co,  X: )  seem 
somewhat  arbitrary.  To  some  extent,  the  value  of  (CO,  A: )  obtained  by  the  CEOF  is  an  energy- 
weighted  average  over  the  spectral  range,  which  accounts  for  their  ratio  approximating  the  given 
c  in  the  above  test  case. 

The  application  of  CEOF  analysis  on  band-passed  time  scries  should  improve  both  the 
reproduction  of  the  time  scries,  by  making  the  signal  more  periodic,  and  the  phase  speed,  by  res¬ 
tricting  the  choice  of  (co,  k )  to  a  narrower  frequency  band.  The  treatment  of  separate  frequency 
bands  is  easily  accomplished  by  constructing  U  (equation  2)  only  over  a  desired  frequency  range. 
In  the  next  example,  the  same  Gaussian  wave  train  is  analyzed  but  in  5  different  frequency  bands 
(0.0  -  0.5,  0.5  -  1.0,  1.0  -  1.5,  1.5  -  2.0,  2.0  -  2.5  cpd).  The  values  of  mode  1  (co,  k)  obtained  for 
each  frequency  band  (Fig.  6)  more  accurately  reproduce  the  given  phase  speed  (represented  by 
the  straight  line).  The  value  of  (co, k)  obtained  from  the  analysis  over  the  entire  frequency  range 
is  shown  for  comparison.  The  time  scries  at  each  station  are  reproduced  from  mode  1  in  each 
band  and  summed  together.  The  reproduced  time  scries  arc  now  more  consistent  with  the  origi¬ 
nal  signal  (Fig.  7),  and  explain  more  of  the  variance  at  each  station  (Fig,  8). 

The  treatment  of  narrower  frequency  bands  dearly  improves  the  analysis,  as  it  must 
since  the  coherence  between  stations  approaches  1.0  as  the  bandwidth  approaches  zero.  For  real 
data,  however,  a  narrower  bandwidth  will  degrade  the  CEOF  results  as  statistical  errors  arise  in 
the  determination  of  the  covariance,  or  equivalently  cross-spectral,  matrix.  To  investigate  this 
effect,  statistical  noise  is  added  to  the  cross-spectral  matrix,  C ,  for  the  five  frequency  bands  dis¬ 
cussed  above  before  applying  the  CEOF.  Following  Davis  and  Rcgicr  (1977),  we  first  add  a 
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random  amount  of  instrument  noise  at  each  station  that  is  not  more  than  10%  of  the  propagating 
variability  and  is  independent  between  stations.  Next,  we  account  for  statistical  sampling  errors 
by  assuming  that  instrument  noise  is  negligible  and  that  the  statistics  of  the  Fourier  coefficients 
used  in  computing  C  have  approximately  joint  normal  probability  distributions.  If  the  process  is 
stationary  and  the  sample  statistics  are  reasonably  accurate,  the  error  matrix,  V,  of  the  cross- 
spectra  is  given  by 

V  =<ClJCnm>^M-lCimCn]  (14) 

where  M  is  the  number  of  independent  spectral  components  used  to  compute  the  observed  cross¬ 
spectra,  and  C  is  the  true  cross-spectra.  In  practice,  V  is  obtained  by  substituting  the  observed 
cross-spectra  for  C  and  M  is  taken  to  be  the  number  of  Fourier  coefficients  in  the  band.  The 
errors  for  C  arc  represented  as  a  linear  combination  of  the  eigenvectors  of  V .  The  coefficients  of 
each  eigenvector  are  obtained  by  selecting  a  random  number  from  a  Gaussian  population  with 
zero  mean  and  variance  given  by  the  corresponding  eigenvalue  (Long  and  Hassclmann,  1979). 
CEOFs  are  then  repeated  for  100  different  noisy  cross-spectral  matrices  with  M  =  18  for  each 
frequency  band.  Phase  speed  is  then  obtained  from  the  five  (eo,  k )  pairs.  The  results  confirm  that 
the  addition  of  instrument  and  statistical  noise  signficantly  alters  the  CEOF  results.  The  average 
phase  speed  is  off  by  7%  with  a  statistical  uncertainty  of  16%.  In  addition,  the  amount  of  vari¬ 
ance  explained  at  each  station  is  degraded  by  up  to  10%  over  the  band-passed  test  case  with  no 
noise,  although  still  better  than  the  original  test  case  over  all  frequencies  (Fig.  9).  The  effect  of 
noise  degradation  (assuming  that  the  process  is  truly  Gaussian)  is  underestimated,  however,  in  the 
specification  of  M .  The  choice  of  M  assumes  a  white  spectra  with  equal  contributions  from  each 
spectral  component  used  to  obtain  the  cross-spectra.  This  is  ah  invalid  assumption  for  the  red 
spectra  of  this  signal;  thus  V  should  be  even  larger.  Prewhitening,  although  not  considered  in 
this  study,  can  make  M  equal  to  the  number  of  spectral  components. 


4.  SUMMARY 

Test  runs  on  synthetic  data  indicate  that  complex  empirical  orthogonal  function  analysis 
is  unable  to  accurately  detect  wide-banded  propagating  signals  similar  to  the  coastal-trapped 
waves  observed  off  the  Pacific  coast  of  Mexico.  In  the  absence  of  noise  and  any  band-pass  filter¬ 
ing,  the  CEOF  recovers  a  propagating  mode  that  is  30%  off  in  phase  speed  and  that  explains  a 
decreasing  amount  of  the  variability  towards  the  ends  of  the  sensor  array.  The  results  arc 
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significantly  improved  by  analyzing  narrower  frequency  ranges.  A  trade-off  arises,  however, 
between  the  CEOF’s  inability  to  handle  wide-banded  signals,  and  statistical  noise  degrading  die 
results  when  the  bands  become  too  narrow.  For  the  ease  of  18  degrees  of  freedom,  which  is  typi¬ 
cal  of  die  degrees  of  freedom  available  in  a  band-by-band  analysis  of  the  actual  observations,  the 
method  returns  a  16%  uncertainty  in  the  phase  speed  and  explains  on  the  order  of  1 0%  less  of  die 
variability  at  certain  stations.  Given  the  shoncomings  of  die  method  illustrated  by  this  very 
idealized  test  ease,  we  abandon  CEOFs  as  a  means  of  detecting  the  observed  coastal-trapped 
wave  signal  off  Mexico. 
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Figure  3. 

The  percentage  of  variance  explained  by  the  CEOF  mode  1  at  each  station. 
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Figure  9. 

The  percentage  of  variance  explained  by  the  combined  mode  1  for  the  five  frequency  bands 
shown  in  Fig.  6  with  statistical  and  instrument  noise  added  to  the  cross-spectral  matrices. 


